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1. Introduction 

SPEC Incorporated was awarded a two-year contract on 5 June 1995 to: 

. Design, build and install a new cloud microphysics instrument on 
the NASA DC-8. 

. Collect data with the new instrument in the SUCCESS field 
experiment planned for the spring of 1996. 

. Develop preliminary software and do quick-look processing of data 
collected by the new instrument. 

On 16 May 1997, the original contract was modified to: 

. Provide additional software to process data from the new 
instrument. 

• Analyze SUCCESS data and publish results. 




Additional milestones were added to the modified contract as shown 
below: 


1) Completion of ice crystal classification of location, concentration, size, 
shape, ice water content and location in cirrus and contrails for 
SUCCESS flights on 4/20/96 and 5/12/96. 

2) Completion of ice crystal classification of location, concentration, size, 
shape, ice water content and location in cirrus and contrails for 
SUCCESS flights on 4/30/96, 5/2/96 and 5/8/96. 

3) Completion of manuscript for submission to Journal of G ®°P h y sic ^| 
Research on the formation of ice crystals and how they are distributed 

in contrails. 

4) Compilation of data base of scattering phase functions and ice crystal 
habits from the Colorado State University cloud chamber. 


5 ) Completion of correlation of individual ice crystal habits scattering 
phase functions from SUCCESS flights on 4/20/96 and 5/12/96. 

6) Completion of correlation of individual ice crystal habits and scattering 
phase functions from SUCCESS flights on 4/30/96, 5/2/96 and 5/8/96. 


7) Completion of review of current parameterizations of radiative transfer 
from cirrus and contrails used in GCMs and development of new 
candidate parameterizations based on analysis of scattering p ase 
function and crystal habit from cloud chamber and SUCCESS data. 


8) Completion of numerical simulations and evaluation of candidate 
parameterizations to be used in GCMs. 

Milestones 1) - 6) were completed and results can be found in 32 monthly 
reports and an interim report submitted on 21 July 1997. The comp e 10 
Milestones 7) & 8) was delayed slightly and the results are included in this fina 
report. In addition, the final report contains a brief summary of work complel ted 
since 1995, a list of publications resulting directly from this rese f r p h a ' S h ° 
papers and conference proceedings that are an extension of this research^ 
Finally Appendix A contains reprints of three journal articles, two of which ^ are 
from the SUCCESS special issue of Geophysical Research Leffa ?’ [L^al 
that is from the NASA FIRE 111 special issue of the Journal of Geo P^ f' 
Research The JGR paper is a report of microphysical properties of Arct^c 
strafus and cirrus clouds' based largely on analysis of data colle^the 
SPEC cloud particle imager (CPI), an outgrowth of the particle imager developed 



for the SUCCESS project. Also, we anticipate submitting the results from 
investigations of Milestones 7) & 8) for publication in a scientific journal. 


2. Summary of Results (Milestones 1-6) 

The new instrument developed for the SUCCESS project was called a n- 
Nephelometer, because it simultaneously performed particle imaging (i.e., Pi or 
the Greek letter n) and measured the light scattered around the imaged particle. 


The initial SPEC proposal which generated the original (two-year) contract 
contained a time frame where the instrument would be developed in the first year 
and flight tested in the second year. However, the SUCCESS field project was 
subsequently scheduled to take place ten months after the contract was signed. 
Both NASA and SPEC scientists acknowledged the desirability of having the n- 
Nephelometer available for the field project, so it was agreed that SPEC would 
attempt to accelerate the development of the new instrument. It was understood 
at the time that the development effort would be streamlined, and that there were 
risks involved, with the possibility that the n-Nephelometer would not be 
available in time for SUCCESS. There were substantial challenges facing both 
NASA and SPEC engineers, including: 


1 The probe and digital data acquisition system were unique. There was no 
existing instrument which produced digital particle images with 5 urn 
resolution, simultaneously measured the scattering phase function and 
recorded all of the data on a high-rate digital data acquisition system. 


2 . 


The probe could not be installed in a standard PMS canister and re 9 u |[® d a 
special mount at the zenith-1 location on the DC-8. NASA and SPEC 
engineers designed a special mount which was fabricated by NASA. Also, 
SPEC had no previous experience with the installation of sensors or 
instrument racks on the DC-8, so new ground needed to be broken in this 


area. 


The overall effort to build and install the n-Nephelometer on the DC-8 and to 
operate it during the field project was successful (Figure 1). However, the 
instrument and data system were literally completed within a few days of the last 
day the DC-8 was available for installing equipment. As a result there were 
aspects of the instrumentation which did not function perfectly. This was 
expected due to the highly accelerated schedule under which the instrument was 

designed and built. 


The optical imaging system, which was new and previously untrie , 
performed exceptionally well and high-resolution (5 pm) digital images of ice 
crystals were recorded. The digital data system was designed to operate at a 
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Figure 1. (top) The n-Nephelometer installed on the NASA DC-8 and (bottom) 
examples of rosettes imaged during the SUCCESS field project, Apnl-May of 1996 








maximum (asynchronous) rate of 25 frames s' 1 , however, duetto a bug in a 
software driver, the maximum obtainable rate was 5 frames s . Thus, wni e 
image quality was excellent, the quantity of image data was less than 
anticipated. However, even though the software bug limited the rate of image 
data the instrument still counted all particles larger than about 20 p.m whil® the 
imaging system was inactive. Thus, it was possible to statistically enhance 
derived quantities such as particle concentration and ice water content. 
Analysis of SUCCESS microphysical data collected in contrails and scattering 
phase functions for particles typically observed in contrails are reported in 
Heymsfield et al. (1998) and Lawson et ai. (1998). Reprints of these articles are 

found in Appendix A. 

The particle imaging system and the digital data acquisition system were 
substantially improved after SUCCESS as a result of parallel funding provided 
by the National Science Foundation (NSF). The improvements were made to an 
instrument built for NSF, called the cloud particle imager (CPI). The CPI was 
installed on the NCAR C-130 for the FIRE III arctic field program conducted May 
- July 1998. Results from the FIRE III field project are reported in Lawson et al. 

(2000), which is included in Appendix A. 

3. Results from Milestones 7 and 8. 


3.1 Introduction 

The impact cirrus clouds and aircraft contrails have on earth’s radiat ion 
balance is poorly understood. Cirrus clouds regularly cover 20% to 30 /o of the 
earth (Warren et al. 1986). Contrails can grow into cirrus clouds, and can 
expand to cover the entire sky if conditions are favorable. Radiative modeling 
has shown that the radiative properties of cirrus cloud particles are extreme 
important (Fu 1996). Detailed measurements of crystal size, shape, and 
radiative properties are necessary to better understand the impact of cirrus 
clouds and contrails (Stephens 1989). 

Model simulations by Fu (1996) showed that optically thin cirrus clouds 
with high concentrations of small particles (less than 10 microns) can have a 
cooling 9 effect at the surface. This is because, visible light energy is reflected 
while infrared energy upwelling from the surface passes through ‘he ctoud since 
small ice crystals act as Rayleigh scatterers for longwave radiation (Wallace and 
Hobbs 1977). This results in a net out-flow of energy. Understanding t 
radiative properties of cirrus particles would improve understanding of the 
radiative impact of cirrus clouds and contrails. 

SPEC designed and fabricated a new instrument, called a Particle 
Imaging n-Nephelometer, and operated it on the NASA DC-8 during the 



Subsonic Aircraft: Contrail and Cloud Effects Special Study (SUCCESS) project. 
The n-Nephelometer imaged cloud particles and made measurements of the 
scattering phase function. Due to instrument problems, the n-Nephelometer 
was unable to take scattering phase function measurements during the field 
project. The scattered light measurement system in the n-Nephelometer was 
repaired and the instrument was later taken to the Colorado State University 
(CSU) cloud chamber. Ice particles similar to those previously imaged in cirrus 
and contrails during the SUCCESS project were generated and scattenng phase 
function measurements were collected (Lawson et al. 1998). 


In Section 3.2 of this report, the scattering phase function measurements 
collected in the CSU cloud chamber are described and analyzed A review of 
common atmospheric model parameterizations is made and new 
parameterizations based on the n-Nephelometer data are deveioped and 
compared It was found that discrepancies between parameterizaters calculated 
from the n-Nephelometer data and commonly used parameterizations were 
greatest for small particles. A cirrus cloud model with the newly created - 
Nephelometer parameterizations is run in Section 3.3 of this report Results 
from parameterizations commonly used in the literature are compared to model 
runs using parameterizations based on n-Nephelometer measurements. 


3.2 Development of New Parameterization (Milestone 7) 


3.2.1 Review of Parameterizations found in the Literature 


Climate and weather 
parameterizations based on the 
parameters commonly used in 
Nephelometer measurements. 


prediction models define clouds with 
cloud’s properties. The following is a list of 
models and how they are related to the fl- 


optical Depth: Optical depth is simply the measure of the cloud s 
optical thickness. The n-Nephelometer was not designed to make 
gross cloud optical depth measurements. A microphysical optical 
depth can be calculated if one knows the concentration and projected 
area of the cloud particles. The n-Nephelometer measures projected 
area exceptionally well and is therefore useful in calculating a 
microphysical optical depth. 


Single Scattering Albedo: The single scattering albedo is the ratio of 
scattered light vs. total extinction. The laser used by the Scatt ® r ® 
Light System (SLS) has a visible wavelength and therefore, very little 
absorption. The single scattering albedo is therefore very close to 1 .0. 



• Asymmetry Parameter: The asymmetry parameter is a measurement 
of how much light is scattered by a particle in the foryvard direction 
versus the backward direction. The asymmetry parameter is 
calculated from the scattering phase function P(0) which was 
measured by the n-Nephelometer. 

Optical Depth 


Microphysical optical depth is shown in eq 1 . 


r = A E * N 


Rg. 1 


T is the optical depth, A E is the extinction per particle (in the visible A E is twice 
the projected area of the particle) and N is the number of particles in the columa 
The more conventional definition of optical depth, based on Beers Law, 

defined as: 


r = — In 


r r 

\hj 


Eg. 2 


where / is the observed radiative intensity and l 0 is the incident intensity^ 
Equation 2 is a bulk optical depth measurement. These definitions converge at 
small optical depths (less than 0.5) where multiple scattering is less importan . 

Asymmetry parameter 

Numerous experimental and theoretical studies have been carried out to 
determine the value of the asymmetry parameter, g, for ice clouds. Table * } 
shows a summary of some of the research efforts as well as results from h s 
work For ice crystals, ray tracing models have been run using basic particle 
shapes (e g Takano and Liou 1995). More recently, complex crystals such as 
aggregated and planar polycrystals, as well as ice crystals with inclusions, have 
also been studied (e.g., Mischenko and Macke 1999). 

The literature contains few experimental measurements of g. jWrerrfland 
satellite measurements have yielded values of 0.7 to 0^8 (Gerber et 
Laboratory measurements have indicated values of g from 0.82 to 0.85 ^ (from 
data in Volkovitsky et al. 1979 and Sassen and Liou 1979). Francis (1995) 
noted that these measured values agree well with simulations, but also noted 
that errors in measuring the forward scattering peak could cause considerable 
uncertainties when calculating g from measured phase functions. 




Ray Tracing results from Takano and Liou 1995 


Particle Shape 

Asymmetry 

Parameter 

Particle Maximum 
Dimension 

Bullet Rosettes 

0.786 and 0.831 

120 and 240 urn 

Equal projected area solid columns 

0.814 and 0.834 

200 urn 

Hollow Columns 

0.843 to 0.858 

L = 300 urn 
hollow 0-75 tarn 

Dendrites 

0.820 to 0.842 

1 0 and 40 pm 

Solid Plates 

0.814 

40 urn 

CaDDed Column 

0.865 

250 pm 

Ray Tracing results from Macke et al. 1993 

Koch-fractal 

0.740 

1 000 urn 

Macke and Mishchenko 1999 

C.rvstal with inclusions 

As low as 0.550 


Parameterizations 

Plates (Mitchell et al. 1996) 

0.85 to 0.95 

30 to 2000 urn 

Columns (Mitchell et al. 1996) 

0.75 to 0.87 

30 to 3000 pm 

Generalized (Fu 1996) 

0.76 to 0.83 

20 to 200 urn 

Aircraft anc 

Satellite measurements ^ 

CIN Probe Gerber et al. (2000) 

0.737 

Glaciated cloud 

Satellite Wielicki et al. 1990 

0.70 

Landsat Re = 60 pm 

I ahoratory measurements/calcula 

tions 

Volkovitsky et al (1979) 

0.85 


Sassen and Liou (1979) 

0.82 






Rosette 

Spheroid 

Fu 

Mitchell 

for 

columns 

Mie 

20 


0.834 


0.847 

0.772 

0.736 

0.867 

40 

0 814 

0.813 

0.800 

0.822 

0.778 

0.759 

0.870 

60 


0.813 I 

0.808 


0.789 

0.773 

0.877 

80 


0.812 

0.812 


0.800 

0.783 

0.879 

100 



0.821 


0.801 

0.790 

0.879 

120 



0.822 


0.809 

0.797 

0.880 

140 



0.829 


0.816 

0.802 

0.881 

160 



0.837 


0.823 

0.807 

0.882 

180 



0.830 0.828 

0.811 

0.884 


Table 1 . Compilation of some measured and theoretical values for the 
asymmetry parameter. 


































3.2.2 Instrumentation 

The n-Nephelometer’s scattered light system (SLS) consists of an array 

of fourteen detectors at different angles an ound the s 03 ''®"" 9 ^ 

measurement angles are: 20, 28, 37, 45 53, , . - ’scattered' from 

152 160 A forward scattering detector catches light scattered rrom 

approximately 2.8 degrees to 9.8 degrees. The SLS laser and opUcs produced a 

(nominaly 1001) polarized beam with a rectangular cross-section. The fonrard 

scattering detector was a circular detector with a rectangular dump spot to block 

the SLS laser. 

The SLS was calibrated using measured water drops and Mie theory 
The forward scattering detector was calibrated for the value of £ • M,e cun* a 
11 decrees. Light scattered off the surface or through the surface of a water 
drop can be highly polarized depending on the angle of incident or refractio . 
^ed" incident onto a water drop it may show much less scattered 
Hoht at certain angles depending on the polarization orientation The in- 

Nephelometer SLS laser was polarized at an angle that could ^make js effert 
sianificant however, the angles where this would cause the most error ' n tn ® 
phase function (60 to 150 degrees) measurement have very little contribu 
(<2%) to the asymmetry parameter calculation. Figure 2 shows a p o o 
theory curve along with the SLS measurement for 20 micron water drops. 

3.2.3 Calculating Asymmetry Parameter from n- 
Nephelometer Data 

the t?rrr irr«^£^ A 

Si"^i a eSmrrj 

with fttetast measured values. Errors introduced by this are d.scussed below. 
The definition of asymmetry parameter is given below. 


n 

g = ^ p (0) cos 6 sin 0d6 la^a 


Where g is the asymmetry parameter, P(9) is the phase 

scattering angle. For our data with 1 degree mcrements equation 3a reduces 

equation 3b. 


180 

g = £p(0)cos0sin0 Eq. 3b 


Before applying eq 3b the phase function needs to be normalized 
according to equation 4: 



Phase Function (arbitrary units) 


Measured and theoretical Phase Function for 



Figure 2 The Phase Function for 20 micron water drops measured 
by the Pl-Nephelometer (gray) compared to the Mie theory calculation 
(black) for 20 micron water drops. 




180 

£/>(0)sin0=l Eq. 4 

After normalization the phase function is multiplied by sin(0)cos(9) and 
then summed to give the asymmetry parameter (Eq. 3b). 


Typical theoretical values of asymmetry parameter range from 0.5 to 0.9 
depending on size and shape of particle. A value of 1.0 would indicate al 
forward scattering, while a value of -1 .0 would be all back-scattering (M.tcheNet 
al 1996). The asymmetry parameter is highly dependant on the forward 
scattering peak. Figure 3 shows a graph of the significance as a function of 
jangle for the phase function for a 20 micron water drop. The significance ^simply 
shows the contribution to the asymmetry parameter from each angular 
measurement as a percentage of total scattered light. 


Mishchenko (personal communication) suggested calculating the 
asymmetry parameter in two parts, the measured phase function part t and the 
diffraction^ part similar to the process described .n Mishchenko and Macke 
(1998) The measured phase function calculation can be treated as the 
Geometric Optics phase function after any diffraction peak is truncated. The 
diffraction component of the asymmetry parameter is the light diffracted by the 
particle in the forward direction. The diffraction asymmetry parameter is always 
very close to 1 .0 in the visible. The measured phase function component of the 
asymmetry parameter is the asymmetry parameter calculated as above except 
the diffraction peak is removed from the phase function. For particles with a aze 
parameter (size parameter = 2rta IX where a is radius and X is wavelength) 
greater than 50 or so, half of the scattered light is from diffraction. The smallest 
parties in this study had a size parameter near 60. Therefore the asymmetry 

parameter is: 


Sd_ Smeas , 

S ““ 2 


Eq. 5 


Where g„ is the diffraction part of the asymmetry parameter and g„«« is he 
asymmetry parameter calculated from the measured phase function (used as the 
geometric optics phase function in Mishchenko and Macke 1998). 

For the n-Nephelometer data the diffraction peak (that came with the Mie 
Theory curve) was truncated at 2.5 degrees. Angles sma Her than . 
were given the same value as the 2.5 degree value. Errors introduced by this 

are discussed below. 



Significance of Phase Function 


Significance of the phase function calculated 
by Mie Theory for 20 micron water drops 



Figure 3. The contribution to the Asymmetry Parameter from each 
angle (Significance = 100*abs[Phase function * sin(theta) * cos(theta))]. 




3.2.4 Asymmetry Parameter Results and Comparison to 
Values Found in the Literature 

Figure 4 shows the asymmetry parameters calculated from the n- 
Nephelometer data. Asymmetry parameter parameterizations calculated from Fu 
(1996) Mitchell et al. (1996) as well as Mie theory are compared to the u- 
Nephe'lometer calculations on Figure 4. The bottom of Table 1 shows the a a 

used to make Figure 4. 

An interesting obsecration is that the smallest particles always have 
asymmetry parameters much larger than the parameterizations^ The smallest 
particles observed by the n-Nephelometer (as well as small particles observed 
in cirrus) were generally spheroidal in appearance. Fu's parameterization is 
based on ray tracing for hexagonal ice crystals. The scattering properties of ic 
spheres and pristine hexagonal plates or columns are quite different (compare 
the Fu and Mitchell parameterizations with Mie theory on fig. 3). The n- 
Nephelometer data shows that small spheroidal ice particles have asymmeby 
parameters closer to those calculated by Mie theory than to he hexagonal 
particle parameterizations Figure 5 shows a typical example °f Part'ce 
observed in cirrus. The particular cloud particles were observed over 0 ^ lah ^ 
by the CPI (the probe developed after the n-Nephelometer). Particle habrt 
classification by concentration as well as radiatively important projected area 
also shown on Figure 5. The small particles are mostly 3 P her °' d f 
appearance. Concentration estimates have shown that cirrus can be made up of 
3m or more particles like these. Contrails observed by the DRI ctaudHopt 
have been shown to be composed of numerous smal P”^"® he ^ ona ' 
crystals. The Fu and Mitchell parameterizations for smaller particles may b 
more appropriate for contrails than for cirrus. 

Possible errors 

The introduction of errors by using Mie theory for 20 micron water 'drops 
to fill in the missing end segments of the measured phase fun®** i were 
thoroughly investigated. Indeed, it seems ludicrous to use theoretical data for 
20 micron spheres to patch in holes for data collected look| h9 ®“P . <- 0 

micron rosettes. To test the reasonableness of mis technique Mie theory 

micron drops (similar projected area to 150 micron bullet rosettes) was 
substituted where Mie theory for 20 micron drops had been used P r e vl ° u ^ y , 
The error was less than 1% when the asymmetry parameter was calculated 
using eq 5, This result becomes obvious when one n0,l0es H , Pf ™ a " 
difference between Mie theory phase functions between zero and 11 ^S" 3 ® 3 ' 3 
the height and width of the diffraction peak. Figure 6 shows the diffraction peak 
for the different Mie theory curves. The diffraction peak for the 50 micron drops 



Asymmetry Parameter 


Asymmetry Parameter calculated from Pl-Nephelometer 
Phase functions compared to theoretical values 


MieTheo 



Spheroid 
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/ 
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20 40 60 80 100 120 140 160 180 

Maximum Particle Dimension foim) 

Figure 4. Asymmetry parameters calculated from the p ^ Ne . phelom [ ?f[ 
data plotted for observed crystals of different sizes and shapes. Data 
are compared to Mie theory, and the Fu and Mitchell parametenzations. 
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Forward Scattering peak for 20 and 50 micron drops (Mie theory) 
and Irregular Ice Particles (Geometric optics) 

1 1 GO Phase function (irregular particle): -* * 

20 micron drops Mie Phase function: 

1000000 \ 50 micron drops Mie Phase function: 


§ 100000 


£ 10000 


1000-1 


1 5 1 2 3 4 5 6 7 8 9 10 11 

Angle (degrees) 

Fimira 6 Mie Theoiv calculated phase function from zero to eleven degree for 
2(?and 50 microifwater drops and the geometric optics (GO) phase function for 
irregular particles from Mishche nko et al. (1996) . 




is much higher at angles less than 1.5 degrees than curve for 20 micron drops. 
The use of the phase function from Mie theory for 20 micron spheres in the 
backscattering region caused almost no change in the calculated asymmetry 
parameter compared to using zeros or interpolating to four times the 160 value 

at 180. 


The degree of truncation of the diffraction peak is a possible source of 
error. To test the magnitude of possible errors the diffraction peak was cut off a 
several different angles from 2.5 degrees to 9.5 degrees. This had little effect 
on the asymmetry parameter (less than 2%). One would expect this to be the 
case from looking at the diffraction peak in Figure 6. The Geometric Optics 
asymmetry parameter for an irregular ice particle (from Mishchenko et aL 19 ^) 
scaled to the other curves is also shown in Figure 6. This confirms that the 
truncation procedure does not cut off any of the Geometric Optics scattered 
light. The Mie theory phase functions slopes are relatively flat in the 2.5 to 11 
degree region similar to the Geometric Optics phase function. 

3.2.5 Summary and Conclusions (Milestone 7) 

The n-Nephelometer cloud probe design by SPEC Inc gathered 
scattering phase function data during laboratory measurements in the CSU 
cloud chamber. These phase functions have been used to calculate the 
asymmetry parameter. A technique was adapted that allows an accurate 
calculation of the asymmetry parameter even when complete information abo 
the forward scattering peak is unknown. Results of this technique applied to th 
n-Nephelometer data show little variation (less than 3%) when the technique is 
purposely applied erroneously. 

The asymmetry parameters calculated from the n-Nephelometer data 
agree reasonably with other laboratory measurements and theoretical 
calculations for larger particle sizes (> 50 pm). For sma\\er particle sizes the 
asvmmetry parameter values were generally higher than current 
parameterizations. This could be due to the fact that the small Particles 
observed were spheroidal rather than the pristine hexagonal shapes which ray 
tracing models are based on. Small particles observed in cirrus clouds are 

generally spheroidal in appearance. 

For accurate calculation of the asymmetry parameter from data with poor 
forward scattering measurements, the asymmetry parameter should be 
calculated in two parts. The diffraction part of the asymmetry P ar ^eter is 
always 1 .0 in the visible when the size parameter is large. The refining part of 
the asymmetry parameter is calculated from the phase function wth any 
diffraction peak truncated. If large areas of the forward peak are not known 
suitable theoretical phase function can be substituted without introducing much 
error (as long as any diffraction-like peak is truncated). The bulk of the error in 



calculating the asymmetry parameter from phase functions is from the height and 
width of the diffraction peak. When the size parameter is greater than 20, the 
extinction efficiency tends towards 2.0. 50% is forward diffraction and 50% is 
either scattered or absorbed. In visible wavelengths where absorption is very 
low, half of the asymmetry parameter is from diffraction and half is from the 
scattering phase function. 


3.3 Numerical Simulations (Milestone 8) 

3.3.1 The model 

The cirrus model used for this study is a Monte Carlo model. The model was 
developed by D. Koracin, V. Isakov and L. Mendez-Nunez (Koracin et al. 
1998). The input parameters for the model include asymmetry parameter, 
single scattering albedo and optical depth for fifteen cloud layers. In addition 
sun angle and surface albedo are also input parameters. The outputs are 
simply radiation flux in the upwards and downwards direction at each level. 

Model Computational Scheme 

Optical properties of the model clouds are characterized by layer optical 
thickness x, single scattering albedo co, and asymmetry parameter g. The 
optical properties of the surface are determined by an albedo A (generally 

assumed to be 18 %). 

For the purpose of this project the model was run with horizontally 
homogeneous cloud layers. For simplicity the photons enter the atmosphere 
vertically When a photon enters the model atmosphere several processes 
are possible. It can be scattered, absorbed, or it can be reflected from the 

surface. 

The photon free-path-length is determined by the optical thickness x: 

r-tw+om 

where k (I) is an absorption coefficient, a{l) is the scattering coefficient, and L 
is the photon trajectory. 

The probability of a photon passing through the atmosphere without 
scattering or absorption is exp(-x). Scattering is considered as a random 
process with a probability ©=a/(K+o). The scattering angle is determine y 



the asymmetry parameter which is used to calculate a probability density of 
scattering angles. 

In order to calculate diffuse radiation fluxes, the Monte Carlo method is 
used by the model. Photons are traced from the top of the model 
atmosphere. Each grid point of the domain has a counter which is 
incremented if a photon passes its grid level in either the upward or 
downward direction. Large numbers of photons (10,000 to 100,000) are 
necessary to achieve sufficient accuracy. 

The model uses the cumulative optical depth as a vertical coordinate. 


f Ik ~ T i-hk + T i,k For ' = 2, n + 1 

= o 

* 

T n+\,k ~ T 0,k 

where / and k represent the vertical and horizontal coordinate in the model 
space. 

For the first trajectory, the counters for all layers are set to zero and the 
photon free path length is determined by: 

a = e r 

r = T-\n(a) 

where a is a random number. Basically, a random number determines how 
much optical depth the photon will pass through before interacting. 

The chance of the photon getting to the surface, being scattered or 
absorbed is then considered. If the free path length enables the photon to 
reach the surface, reflection from the surface is also considered. 

If the photon free path length is less than the total optical depth in the 
simulated cloud, a photon scattering event is considered at the cloud 

position: 

Absorption was not considered for this project (the single scattering albedo 
was set to 1.0 which is reasonable for the visible). A random number 
determines the scattering angle of the photon based on the local asymmetry 
parameter. After a photon is scattered, a new photon free path length is 
determined using the random number generator. Once the photon reaches 



the statistical noise level (as determined by d discussed earlier), it is no 
longer considered. 


Counters are incremented based on how far the photon makes it: 


pr ^ - A . 

r n+ 1 _ A r n+\,k 

p total i _ p total i 
1 i,k i,k 

r? totals rr total T , 17 T 

+ r i,* 


A new photon enters the system once the first is finished and the same 
process is repeated. Fluxes are calculated based on the numbers counted at 

each layer. 


Emphasis for Modeling Runs 


The main difference between the n-Nephelometer calculated asymmetry 
and the common parameterizations (Mitchell et al. 1996 and Fu 1996) are at 
the small particle end. Figure 5 shows that typical cirrus contains high 
number of small spheroidal particles which make up a significant portion of 
cirrus cloud extinction (i.e. projected area). For smaN part'des the 
parameterizations show low values of asymmetry parameter ( . o . ). 

The n-Nephelometer data gave values for asymmetry parameter that were 
much higher for 20 micron spheroidal particles (0.85). The discrepancy 
between the parameterizations and our calculations decreased for larger 
particles. Since the radiatively important small particles showed the most 
discrepancy, the modeling effort was concentrated there. 


3.3.2 Model Results 

Vogelmann and Ackerman (1995) used observed surface radiation fluxes 
to relate to cirrus cloud properties. Cirrus clouds were modeled and the effect of 
small perturbations in the cloud properties was investigated. A similar technique 
is used here to analyze our model results. 

For the first model runs, we fixed the surface flux and varied the cloud 
properties based on Part l results. A cloud with five layers of optical depth! 15 
was modeled. The asymmetry parameter for each of these layers was defined to 
be 0.75 to represent the currently parameterized asymmetry parameter for small 
cirrus particles. 0.75 is about what the Fu (1996) and Mitchell et al. (1996) 
parameterization schemes give for small particles commonly observed in cirrus 
The model was run and the flux at the surface was determined. The asymmetry 
parameter was then changed to 0.85 which is the value calculated frorr » the > n- 
Nephelometer data. A number of model runs were then completed with differen 
optical depths in an effort to match the surface flux. A perturbation in the optical 



depth of about 2 percent was necessary to obtain equal flux at the surface. A 
cloud with an asymmetry parameter of 0.75 produced the same surface flux as a 
cloud with an asymmetry parameter of 0.85 and a 2% different optical depth. 
Table 2 shows the parameters input into the model to achieve equal surface 

flux. 


Asymmetry Parameter 

0.75 

0.85 

Optical Depth (5 layers) 

0.50 

0.4905 

Surface Flux W/m 2 

392 

392 


Table 2. Input parameters for model runs to achieve equal surface fluxes. 
Optical depth and asymmetry parameter were varied. 

The second model test used assumptions based on observed cirrus 
particles and the particles used for asymmetry parameter parameterization^ 
Naturally occurring small cirrus particles are spheroidal in appearance (fig. 4) 
while the particles parameterized by Fu (1996) and Mitchell et al. (1996) were 
hexagonal. The microphysical optical depth can be calculated using the 
average projected area of the particles. Projected area for a sphere and for a 
hexagonal particle (equal to its surface area / 4 for a randomly oriented convex 
particle) were used to determine microphysical optical depths for equal Ice 
Water Content (IWC) clouds. The microphysical optical depth is theoretically 
equal to the cloud bulk optical depth when the optical depths are low. Model 
runs were made on low optical depth clouds for equal IWC clouds with 
appropriate asymmetry parameters and optical depths. Table 3 shows the clou 
particle properties input into the model and the resulting surface ux. 
difference of 19 Watts/m 2 between the clouds was determined by the model. 


Particle shape 

Spheroid 

Hex-Crystal 

Optical Depth 

0.05 

0.0423 

Surface Area pm 2 

314 

371 

Volume (mass)pm 3 

4188 

4188 

Asymmetry Parameter 

0.85 

0.75 

Surface Flux W/m z * 

1069 

1050 


Table 3 Cloud parameters for second set of model runs and the resulting 
surface flux calculation. The surface flux is based on the fraction of solar 
irradiance that falls between 0.25 and 1.41 pm (from Fu 1996). 


3.3.3 Summary and Conclusions (Milestone 8) 

Model results have shown that differences in the asymmetry parameter 
can cause significant differences in optical depth and surface flux. A cloud 





model was run in an effort to find the difference in optical depth that could be 
caused by using an incorrect asymmetry parameter. A cloud with an asymmetry 
parameter of 0.75 was modeled, then the asymmetry parameter was changed to 
0.85 and the optical depth was varied in an effort to match the surface radiation 
flux of the first model run. The optical depth had to be changed by 2% in order 
to get the same surface flux in the second model run. Equal surface fluxes can 
be measured by clouds with two percent differences in their optical depth 
depending on what asymmetry parameter is used. 

A second modeling effort was done to try to compare identical ice water 
content clouds. The cloud particles were designed in an effort to try to match 
observed clouds (spheroidal particles with asymmetry parameter 0.85) to those 
that were used for asymmetry parameter parameterization calculations 
(hexagonal particles with asymmetry parameter 0.75). A microphysical optical 
depth was calculated by estimating twice the projected area for a column number 
of particles. Optical depths were kept low in order to reduce the possibility of 
multiple scattering. A model run was done on each of the equal ice water 
content clouds. The surface flux difference was 19 W/m 2 . This huge difference 
could mainly be accounted for by the different optical depths for the clouds. This 
result also shows that making cloud mass estimates from surface flux 
measurements is highly dependant on assumed particle microphysics. This 1 9 
W/m 2 value should be thought of as an upper limit as it assumed perfect 
hexagonal crystals and perfect spheres as particle shapes. 
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Growth of ice crystals in a precipitating contrail 
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Abstract. . 

This study examines how jet aircraft contrails develop 

precipitation trails, using data collected on 1 ay, 
1996 during SUCCESS. The DC-8 sampled the pre- 
contrail conditions, produced a contrail largely m clear 
air at -52° C, and sampled the contrail and developing 
trails for over an hour. The environment was highly 
ice-supersaturated, reaching nearly water saturation m 
some locations. Inside the contrail core, almost all ice 
particles remained small (~ 1 to 10 fJ.ni) due to lug l 
crystal concentrations (~ 10 1 to 10 2 cm' 3 ) which re- 
duced the vapor density to saturation. Mixing of moist 
environmental air and vapor-depleted contrail air pro- 
duced localized regions of supersaturation along the 
contrail periphery, where crystals grew to several hun- 
dred microns at about 0.1 /xm s . These particles could 
then fall from the contrail into the vapor-rich, unde- 
pleted, supersaturated environment below. As heavier 
crystals left the contrail, others moved into the regions 
of ice supersaturation. Precipitation trails developed as 
this process continued over time. 


2. Data Overview 

The DC-8 made five sampling runs from ~ 10.5 to 
11.5 km (MSL pressure altitude), between 22:38 and 
24 00 UTC Figure 1 identifies these flight tracks as pen- 
etrations one through five (Pens. 1 - 5). A side effect of 
multiple sampling runs was the generation of additional 
contrails. Lidar data, along with the onboard observer 
notes, were used to differentiate ambient cirrus from the 

contrails, when possible. or 

Pen. 1 occurred at Z = 10.4 km and T = -47 C, see 

Fig 1A The aircraft carried a diode laser hygrometer 

(DLH), see Vay et al. (1997), to measure the water va- 
nor mixing ratio which was converted to RH (tig. ) 
based on the air temperature. The RH at ice saturation, 
RHj, was derived from T. 

The DC-8’s rear-facing video camera recorded con- 
trail formation in the wake of the aircraft. A contrail 
was produced during Pen. 2, at 11.1 kin where T 
-52°C see Fig. IB, and RH >> RHj, see Fig. 1C. This 
contrail persisted for more than six hours during which 
v/iciKlp in the (TOES satellite imagery. 


timo it WPS 


1. Introduction 

Condensation trails or contrails are cloud-like stream- 
ers that frequently form behind aircraft hying in co , 
clear humid air. Critical environmental factors for con- 
trail generation are the temperature pressure, and wa- 
ter vapor density, see Schmidt (1941), Appleman (1953) 
or Khrgian and Mazin (1955). Persistent contrails oc- 
cur when the ambient relative humidity with respect to 
water (RH) is at least at ice saturation (RHj). Contrails 
occasionally produce precipitation trails. 

There are few in-situ or remote measurements as- 
sociated with the precipitation trails of contrafls (We- 
ickmann,1949; Knollenberg,1972; Konrad and Howard, 
1974- Schumann, 1994). This paper examines the envi- 
ronmental factors that enable contrail particles to grow 
into larger precipitating particles, using data from 
Mav 1996 when the NASA DC-8 research aircraft gen- 
erated a contrail northwest of San Francisco C A as part 
of the Subsonic aircraft: Contrail and Cloud Effec s 
Special Study (SUCCESS) field program. 
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7 iR i DC-8 measurements on 12 May 1996. Top: DC-8 trade. Pen- 
18-2250 UTC (solid line); 2: 2253- 2313 (dotted); 3: 2315-2332 (long A*as , 

2333 2349 (dark solid and small dots); 5: 2350-2400 (dash dotted). B: Time 

Tof B priure height (km) (left ordinate) and corre.pond.ng temper^e 

ght), with penetration numbers shown, and C. Time plo ( 

(dashed), and locations where NO > 100 ppmv (dots). 
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23:23:00 


23:39:20 23:39:30 23:39:40 23:39:50 23:40:00 

Time (hh:mm:ss) 

froMnni ° nTf"? f ', 0 "‘ Pe "' 3 (“!>)• »>'<■" «'<■ dc-s m 

•4 kin, and Pen. 4, when it descended from 11.5 to 11.0 km. 


23:40:10 


The DC-8 ’s lidar, see Uthe et al. (1998), was nadir- 
pomtmg for the first seven min. of Pen. 2 and detected 
patchy cirrus extending to between 0.5 and 1.0 km be- 
low tin' aircraft. For the remainder of Pen. 2 visual 
observations indicated that the contrail was often dis- 
tinct, but occasionally mixed with the ambient cirrus. 

iiring Pen. 3, the DC-8 entered the Pen. 2 contrail 
at least, sixteen times. Inside the contrail, NO concen- 
trations significantly exceeded the background levels of 
css than about 50 pptv, see Weinheimer et al (1997)- 
we define periods where NO > 100 pptv as “high NO” 
periods to identify the contrail core. Fig. ID shows that 
w heic NO > 100 pptv, RH * RH, = 62%, and where 

■ <100 pptv, outside the contrail, primarily RH >> 

n H l . 


Lidar observations from an ascent through the con- 
trail during Pen. 3 show the contrail as an area of rel- 
atively high brightness (reflectivity) with a cellular-like 
structure that extends through a several hundred meter 
deep layer (Fig. 2A). As the aircraft climbed through 
the contrail layer, the lidar (downward viewing) became 
more visible. The lidar observations show the presence 
of less reflective, non-contrail cirrus below. The sharp 
contrast between contrail and clear air at the base of 
the contrail indicates that virga had not yet developed. 

tiring Pen. 4, the DC-8 penetrated visible contrail 
at least six times, see Fig. IB. Again, RH « RH, inside 
the contrail while RH » RH t outside the contrail. 
Virga from the contrail had begun to develop as noted 
near the base of the contrail. The presence of contrail 
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precipitation was confirmed by the visual observations 

of the pilots. _ 

During the Pen. 5, the DC-8 climbed from 10.5 to 
11.5 km and entered the contrail only during one period. 

RH >> RHn see Fig. 1C. 

3. Microphysical Measurements 

Here we discuss the microphysical data from 233920 
and 224010 when the DC-8 descended at a constant rate 
through the contrail and associated virga; microphysi- 
cal data for other periods is presented in Lawson et al. 
119971 Figure 2B shows the lidar data. Air parce 
trajectory calculations using the DC-8 measured hori- 
zontal winds indicate that this part of the i contrail was 
produced at approximately 230400 and had not been 
repenetrated in the 35 min. from the time of produc- 

The instruments used to characterize the contrail mi 

crophysics include the following: , sp > 

( 11 A Multiangle Aerosol Spectrometer Probe (MAS ) , 
see Baumgardner et al. 1995, measuring size spectra 
from 0.4 to 40 pm diameter; concentrations between > 

1 and > 40 pm are used here; 

(2) A particle imaging nephelometer (PI), see Law - 
son et al. 1997, providing ice particle images and size 
spectra for particles exceeding ~ 20 pm diameter; 

(3) A video ice particle sampler (VIPS), see Heyms- 

field and McFarquhar (1996) to estimate the particle 
size spectra for particles in this study above 30 pm di- 
ameter; and, . , 

(4) A PMS 2D-C imaging probe with a lower detec- 
tion threshold of between 50 and 100 pm. 

Figure 3 shows data taken through the contrail core 
as a function of height during the descent along the up- 
per abscissa and time along the lower abscissa. The air- 
craft descended through the contrail and precipitating 
region at an approximately constant rate so time and 
height are nearly linearly related; temperatures warmed 
from -56 to -52°C. The shaded region in Fig. 3 dehnes 
the contrail core (from NO). Here, RH « RH 4 = 62 to 
64 %, while on either side of the core, RH >> RH,- 
The contrail core contained high concentrations o 
small ice crystals see MASP data, Fig. 3B. The ad- 
jacent, non-contrail regions contained orders of magni- 
tude fewer crystals predominantly < 100 pm. W ie 
maximum particle dimensions were at most <100 pm 
in the contrail core, see Fig. SC, they extended to 250 
pm in the high RH, peripheral regions see Fig. 3C and 
Lawson et al. 1997. While the particles in the con- 
trail core were small and indistinct in shape (see im- 
ages), particles in the high RH regions were typically 
well-defined bullet rosette crystals, typifying growth in 
a high RH environment. 

Note that large crystals are also observed above the 
contrail core in Fig. 3. Gierens (1996) shows that the 
rms air vertical velocity in the contrails are on the order 
of the crystal terminal velocity, suggesting that vertica 
motions could explain transport of crystals above the 
contrail. 


4. Discussion 

The environmental and contrail conditions discussed 
below probably occur in most situations where precipi- 
tation trails are produced by “cold” contrails, i.e. those 
forming below -50°C, like the one studied here. 

. The ice particles in the contrail core primarily develop 
in an environment where RH is close to but slightly 
greater than RH t . As the small particles leave the con- 
trail interior, the lack of natural cirrus is necessary to 
provide relatively high RH environmental air where the 
contrail crystals can grow. 

. A sustained region of high RH air is required because 
the ice crystals grow' slowly at temperatures lower than 
-50 °C see Gierens 1996. The larger crystals in the 
contrail periphery (Fig. 3) grew to at most 250 pm from 



23:39:20 23:39:45 

Time (UTC) 

Fig. 3. Relative Humidity and microphysical properties fo (J3:39.20 to 
■23-40 10 during DC-8 descent, for same time period as in Fig. 2B. Shading 
high NO’ region. Water vapor and microphysical properties as a function of Urn 
during descent through contrai. for time period corresponding to hdar data n 
Fig 2 A: RH, examples of observed ice crystals. B: Conceit from MAS ( 
to 40 pm), PI nephelometer, VIPS, and 2D-C. C: Size of largest particle from 

2D-C. 
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the time of generation approximately 2100 sec earlier, 
or with an average growth rate for the largest crystals 
of 0.12 /im s -1 . 

• A relatively deep layer of high RH air is required for 
ice particles outside the contrail perimeter to prevent 
sublimation and promote further development follow- 
ing ice crystal fallout into the environment below. In 
effect, this deep layer of high RH air makes the precip- 
itation trails visible. Our observations show that the 
environmental conditions were supersaturated with re- 
spect to ice for depths of at least 0.7 km below the 
initial contrail height for most of the contrail, see Fig. 
2. These broad, sustained regions of high RH suggest 
the presence of large scale upward vertical motion. 

Finally, high initial concentrations of ice crystals within 
the contrail core are required because the existence of 
a precipitation trail implies the continuous flux of ice 
crystals from the contrail to the precipitation trail over 
time periods exceeding 30 min. Furthermore, high ini- 
tial ice crystal concentrations are also needed because 
the contrail spreads substantially in the horizontal see 
Gierens 1996. 


plying sustained upward vertical motions), and a deep 
layer of high relative humidity. 
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5. Conclusion 

Occasionally, contrails generate precipitation trails 
such as those observed during the 12 May 1996 SUC- 
CESS contrail study. The unique aspects of this study 
included: 

(1) The generation of a distinctive-shaped contrail 
by a research aircraft equipped with a State-of-the-art 
instrument, complement in a largely cirrus-free environ- 
ment; 

(2) Environmental conditions that were highly super- 
saturated with respect to ice; 

(3) In-situ sampling and lidar mapping of the resul- 
tant contrail over a one hour period; and, 

(4) A distinguishable contrail that persisted for over 
six hours. 

Inside the contrail core, ice crystal concentrations 
were ~ 10* to 10 2 cm -3 , with few crystals > 100 pm 
in dimension. In the peripheral regions, concentrations 
of crystals < 100 pm were one to two orders of magni- 
tude lower, but concentrations of crystals > 50 pm were 
higher, 0.01 to 0.1 cm' 3 , comparable to those measured 
by Knollenberg. Significant concentrations of crystals 
> 100 pm were noted on the contrail periphery, and 
•sizes up to 250 pm were detected. We suggest that 
turbulence allows some ice crystals to be selected for 
preferential growth, serving as a source of particles for 
the associated precipitation trails. 

Particular environmental conditions are shown to be 
needed for contrails to produce precipitation trails, es- 
pecially at the temperatures studied here as ice crystal 
growth rates are very low at temperatures < -50°C. 

In the highly ice-supersaturated growth environment 
studied here, linear growth rates to 250 pm averaged 
only 0.12 pins -1 . Favorable environmental conditions 
include: a largely cirrus-free environment, a sustained 
growth period in high relative humidity conditions (im- 
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Abstract A persistent contrail in the shape of a racetrack was 
generated by the NASA DC-8 research aircraft during the SUC- 
CESS project The contrail was visible on GOES imagery for six 
hours. Microphysical measurements collected by the DC-8 show 
that after 40 min the core of the contrail contained mostly small 
(1 to 20 pm) ice particles in concentrations >1000 L* 1 , with larger 
(>300 pm) ice crystals in concentrations < 10 * LA In contrast to 
the core, the contrail periphery contained about an order of mag- 
nitude less ice particles in the 1 20 pm range and about three 

orders of magnitude more ice particles >300 pm. The larger ice 
crystals in the periphery were mostly columns and bullet rosettes, 
similar to habits of larger ice crystals found in ambient cirrus in 
the area Measurements suggest that the shape of phase functions 
of randomly-oriented columns and rosettes are mostly featureless. 
The measured phase functions are closest in shape to those pre- 
dicted by ray-tracing theory for random-fractal and spatial- 
dendrite ice crystals. 


1. Introduction 

A persistent contrail was formed by the NASA DC-8 research 
aircraft at about 1 1 km (-55°C) offshore of San Francisco on 12 
May 1996 during the Subsonic Aircraft: Contrail and Cloud Ef- 
fects Special Study (SUCCESS) project Heymsfield et at (1997) 
(in this issue) show the racetrack-shaped flight track of the DC-8 
and the resulting contrail, along with measurements of lidar back- 
scatter, relative humidity, particle concentration and size. The 
persistent contrail, with its distinctive racetrack pattern, was visi- 
ble on GOES satellite imagery for about six hours while it moved 
eastward over the California mainland and dissipated in the lee of 
the Sierra Nevada mountains. Precipitation streamers, visible 
from the bottom of the contrail, were occasionally noted by the 
pilots on the DC-8. 

In situ measurements from Heymsfield et al. (1997) show that, 
in the ambient environment and at the periphery of the contrail, 
the vapor content of the air was typically 20 - 40% supersaturated 
with respect to ice. On the other hand, the core of the contrail 
was at vapor saturation, the excess vapor being depleted by small 
(< ~20 pm) particles in concentrations of about 1 to 10 cm 3 . The 
measurements suggest that a turbulent process results in a small 
percentage of the particles being mixed from the core into the va- 
por-rich environment at the periphery of the contrail, where they 
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grow to sizes > 300 pm and may form virga (i.e., precipitation 
streamers). The in situ measurements are in general agreement 
with earlier work by Knollenberg (1972) and Heymsfield (1973). 

Persistent contrails have been observed to coalesce into clouds 
which take on the appearance of cirrus and can act to measurably 
reduce the incoming solar radiation in heavily trafficked flight 
corridors of Europe and elsewhere (Bakan et al. 1994). The light- 
scattering properties of ice crystals and radiative characteristics of 
cirrus clouds are strongly dependent on particle shape (Mish- 
chenko et al. 1996a). A new in situ sensor, which was installed 
for the first time on the DC-8, provided high-resolution (5 pm) 
digital images of the ice particles. The instrument, called a parti- 
cle imaging nephelometer, or n-Nephelometer, simultaneously 
records the particle image and measures the scattering phase 
function of the imaged particle. The shapes of ice particles wdth 
images >—50 pm are readily discernible and provide information 
of typical particle shapes in the periphery of the contrail and in 
cirrus which occurred naturally in the vicinity of the contrail. 
Measurements collected in the Colorado State University (CSU) 
dynamic cloud chamber (DeMott et al. 1 990) after the SUCCESS 
project were used to correlate the imaged particles with the scat- 
tering phase function. 

The measurements on 12 May 1996 show that the size distri- 
butions of particles within the core of the contrail differ from 
those on the periphery, in the precipitation streamers and in the 
ambient cirrus. Consequently, the light-scattering properties of 
these regions can be expected to differ, which in turn influences 
their optical and radiative characteristics. 

2. Instrumentation 

The NASA DC-8 was extensively instrumented for micro- 
physical research. In addition to the n-Nephelometer, micro- 
physical measurements from instruments discussed in this paper 
include: 

(1) A multi-angle aerosol spectrometer probe (MASP), de- 
scribed by Baumgardner et al. (1996), measures particle size dis- 
tribution in the range from 0.3-40 pm; 

(2) a video ice particle sampler (VIPS), as described by 
Heymsfield and McFarquhar (1996), uses 30 Hz video images to 
estimate the particle size spectra for particles of at least 5 pm di- 
ameter. During this flight, there was a focusing problem which 
limited the minimum detectable size to 30 pm; and 

(3) a PMS 2D-C probe (Knollenberg 1981) with a lower de- 
tection threshold of between 50 and 100 pm. 

The n-Nephelometer casts an image of a particle on a 400,000 
pixel solid-state camera by freezing the motion of the particle 
using a 25 ns pulsed, high-power (60 W) laser diode. Unique op- 
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tical imaging and particle detection systems (PDS) precisely de- 
tect particles and define the depth-of-field (DOF) so that at least 
one particle in the image is almost always in focus. This elimi- 
nates out-of-focus sizing errors that have plagued the conven- 
tional 2D imaging probes (Korolev et al. 1991, 1998). 

The n-Nephelometer was designed to record up to 25 frames 
s*\ although a software bug prevented the instrument from re- 
cording more than 5 frames s l . This severely limited the imag- 
ing sample volume. 1 However, the particle detection system was 
always active and recorded the number of particles > - 20 pm 
which traversed the 3 mm by 4 mm particle detection area. At 
225 m s\ the sample volume of the PDS is 2.8 L s* 1 . Particle size 
distributions were computed by averaging several seconds of im- 
age data and scaling the concentration using measurements from 
the PDS. 

The scattered light system (SLS) consists of twenty-eight 1- 
mm optical fibers connected to microlenses bonded on the surface 
of avalanche photo diodes (APDs). The fibers are placed at dis- 
crete angles of 11, 20, 28, 37,45, 53,61,70, 110, 118, 127, 135, 
143, 152 and 160°. The field of view of each fiber is about ± 2°, 
except for the 1 1 0 fiber which integrates forward scattered light 
from about 3.5 -> 1 8°. 

3. Size distributions and shapes of ice crystals 

Heymsfield et al. (1998) show five DC-8 racetrack-shaped 
flight patterns on 12 May 1996 which often overlap and some- 
times intersect ambient cirrus clouds which were scattered over 
the area. The ambient cirrus and overlapping contrails con- 
founded unique determination of the origin of the ice particles. 
The central cores of the contrails were often detectable from dis- 
tinctly elevated NO* concentrations (Weinheimer et al. 1997). 
The peripheries of the contrails were delineated on either side of 
the central core. The contrails and precipitation streamers were 
occasionally visible from the cockpit of the DC-8, and intercom 
voice notes were used to determine the position of the aircraft 
relative to the visible contrail. Lidar observations, shown in 
Heymsfield et al. (1998), were also used to identify cirrus that 
was above or below the aircraft. These composite observations 
were used to select representative regions of the contrail core, pe- 
riphery and ambient cirrus. 

Ambient cirrus. The cirrus clouds in the area near the con- 
trail were mostly thin, patchy and contained occasional particles 
up to 1 mm in maximum dimension. Fig. 1 shows a combined 
particle size distribution using the MASP, VIPS, ri-Nephelometer 
and 2D-C probes along with examples of images of particles ob- 
served by the ri-Nephelometer from 225230 - 225630 GMT. 
There is good general agreement in the shapes of the particle size 
distributions measured by the VIPS, IT-Nephelometer and 2D-C 
probes for particles from about 75 - 500 pm. In Fig. 1, the col- 
umns extend to about 200 pm in length and the maximum dimen- 
sions of bullet rosettes are up to 600 pm. The shape of the small- 
est (< 50 pm) particles often cannot be determined They usually 
appear to be mostly spheroidal in shape, however, this may be 
due to the II-Nephelometer 5 pm pixel resolution, which is not 
always adequate to resolve the edges of plates < 50 pm. As a re- 
sult, in this paper we have elected to classify small particles that 
appear spherical in shape as “unresolved spherical”. The shape 
of particles from 50 to 200 pm are mostly columnar and/or ro- 
settes, and particles > -200 pm are usually bullet rosettes. 


' Upgrades after completion of the SUCCESS project included a 1 mil- 
lion pixel CCD camera, 2.3 pm pixel resolution, increased optical contrast 
and 40 frames s ' data recording rate (see Lawson et al. 1 998). 


^ 1E+4 

f IE+2 

a lE+0 

I 

I ,« 



0 100 200 300 400 500 600 

Diameter (pm) 

Figure 1 . Particle size distributions from the MASP, VIPS, 
IT-Nephelometer (Pi), 2D-C measurements and composite distri- 
bution in (a) ambient cirrus, (b) contrail core and (c) contrail pe- 
riphery. Examples of images from the n-Nephelometer are 
shown next to each case. 


Contrail cores. The central cores of contrails typically con- 
tain high concentrations of small ice particles (Knollenberg 1972; 
Heymsfield et al. 1998). A contrail core, identified by NO, > 100 
ppt was sampled from 233938 to 233952 (Heymsfield et al. 
1998). The composite size distribution shown in Fig. 1 is domi- 
nated by the small particles seen by the MASP. Images of the 
small particles recorded by the n-Nephelometer were classified 
as “unresolved spherical”. High-resolution replicas seen on the 
Desert Research Institute (DRI) home page (http://www.dri.edu) 
suggest that particles < 20 pm in contrails are often plates and 
columns, combined with irregular and spheroidal shapes. How- 
ever, the crystals seen on the DRI home page were all observed in 
an environment with low (< ice saturation) relative humidity and 
may not be representative of the 12 May contrail, which was typi- 
fied by very high relative humidity. More than 99% of the ice 
particles in the core of the 12 May contrail were < 20 pm and the 
low concentration of 100 to 250 pm rosettes may have mixed in 
from the contrail periphery (Heymsfield et al. 1998) 

Contrail periphery and precipitation streamers. 

Heymsfield et al. (1998) show data from a penetration of a 
contrail that was about 40 min old. The DC-8 penetrated the 
contrail downward from above. Data from the lidar suggested 
that the contrail itself was not contaminated from cirrus aloft and 
that there were some regions that contained precipitation stream- 
ers. Heymsfield et al. (1998) show MASP measurements of total 
concentration of particles from about 1 to 10 pm in the contrail 
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periphery that are 1 to 2 orders of magnitude lower than in the 
core. The composite particle size distribution and II- 
Nephelometer images shown in Fig. 1 suggest that the contrail 
periphery contains particle types that are very similar to ambient 
cirrus, with the exceptions that there are additional small particles 
and fewer particles in the 300 — 600 pm range. Heymsfield et al. 
(1998) suggest that a turbulent process occurs whereby some par- 
ticles in the core are mixed into the vapor-rich region of the pe- 
riphery and grow to larger sizes. 

As the DC-8 continued its descent, it may have also penetrated 
precipitation streamers below the contrail*, however, some of the 
regions also contained ambient cirrus and separating the two was 
unambiguous. Basically, the data suggest that the particles con- 
tinued to grow in the vapor-rich environment as they precipitated 
from the contrail, sometimes reaching sizes up to about 1 mm. 
The habits of particles in the precipitation streamers appeared to 
be dominated by bullet rosettes. 


4. Light scattering properties of the ice crystals 

Calculations shown in Mishchenko et al. (1996a) suggest that 
using the wrong scattering phase function in retrieving cloud op- 
tical thickness can result in an overestimation or underestimation 
of optical thickness by more than a factor of three. The 11- 
Nephelometer was used to measure the phase functions of ice 
crystals in the CSU cloud chamber. Before the data collection, 
extensive laboratory tests were conducted with water drops and 
glass beads. Comparisons with Mie scattering theory for 20 to 
200 pm water drops and IT-Nephelometer measurements were in 
good relative agreement. However, measurements of 50 and 100 
jjjn glass beads scattered nearly an order of magnitude more light 
in the region from 80° to 120° than predicted by Mie theory for a 
spherical glass bead Gayet (1997 - personal communication) ob- 
served a nearly identical departure from Mie theory in the 80° to 
120° region when measuring phase functions from glass beads 
with a polar nephelometer (Gayet et al. 1997). Mishchenko et al. 
(1996b) found very similar results when comparing micron-size 
spheres with equivalent area spheroids and soil particles. This 
suggests that even minute degrees of non-sphericity can lead to 
large differences in phase functions. 

Fig. 2 shows phase functions for 20 pm water drops generated 
in the laboratory and ice particles observed in the CSU cloud 
chamber along with examples of n-Nephelometer images. All of 
the phase functions shown in Fig. 2 are normalized to the light 
intensity measured at the 20° scattering angle for a 20 pm water 
drop. The dynamic chamber was operated in a way such that a 
cloud was formed by an adiabatic expansion at about -20°C (540 
mb), seeded at -28°C and ice crystals grew as the expansion con- 
tinued to about -42°C (380 mb). The 11-Nephelometer images 
were sorted by size and the particle shapes observed in the 12 
May contrail and ambient cirrus. Figure 2 shows that the relative 
intensities of the phase functions for small (20 to 50 pm) particles 
are strongly a function of particle shape. The data show that in 
the 20 to 50 pm (Sm.) size range, progressively more light is 
scattered (from about 40° -► 160°) as the particle shape becomes 
more complex, from water drops — > unresolved spheroids -» ir- 
regulars — ► columns and rosettes. In Fig. 2, for the 50 to 90 pm 
(Med.) and 90 to 210 pm (Lg.) size ranges, it appears that the in- 
fluence of particle shape on phase function decreases. For parti- 
cles > 50 pm, the data suggest that size is the significant factor in- 
fluencing the amount of light scattered 

Figure 2 also shows phase functions derived from ray-tracing 
for a dendrite (Takano and Liou 1995), a column and a fractal 
(Mishchenko et al. 1996a). The prominent 22° halo seen for the 



Figure 2. Normalized phase functions derived from meas- 
urements of ice crystals in the CSU cloud chamber compared 
with theory. Lighter lines offset by one decade toward the top of 
the figure are phase functions derived from ray-tracing for ran- 
dom-fractal and columnar crystals (Mishchenko et al. 1 996) and a 
spatial dendrite (Takano and Liou 1 995). Shown in lower part of 
figure (top to bottom) measurements from columns (Col), bullet 
rosettes (Ros), irregular shaped crystals (Irr), and unresolved 
spherical (Unr), 20 pm drops, and Mie Theory (Thry) for 20 pm 
drop. Size bins indicated in legend for measured crystals include 
20-50 pm (Sm.), 50-90 pm (Med.), and 90-210 pm (Lg.). Num- 
ber (N) of measurements for each crystal type is included in leg- 
end. Also, examples of Pi-Nephelometer images of bullet ro- 
settes (top row), columns and irregular shapes (middle row), and 
unresolved spherical crystals (bottom row) are shown in the mid- 
dle right of the figure. 


column is not observed in the measurements. Mostly, the meas- 
ured phase functions are featureless with some suggestion of a 
peak near 145°. The best agreement between the measurements 
and ray-tracing theory appears to be with the fractal and dendntic 
crystals. Nikiforova et al. (1978) and Volkovitskiy et al. (1980) 
also made measurements in a cloud chamber. They used a 
“Rassvet” device (Nikiforova et al. 1978), which averaged meas- 
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urements over a range of discrete angles and a time period of 10 
s. The phase functions shown by Nikiforova et al. (1978) do not 
display the 22° halo for the small ice crystals and the authors 
comment that the 22° halo was only observed in clouds with low 
optical density. Interestingly, Nikiforova et al. ( 1 978) do show a 
slight bump at 145° for both ice and water clouds. Volkovitskiy et 
al. (1980) show a 22° halo in their measurements, although they 
comment that it is not nearly as prominent as predicted from ray- 
tracing. The principal difference between the measurements us- 
ing the Rassvet device and the Il-Nephelometer, is that crystals 
are sucked through a 25 mm sample tube in a random orientation 
in the FI-Nephelometer, and they are allowed to fall undisturbed 
through the Rassvet device. Apparently, under conditions where 
the optical density is low and hexagonal crystals are allowed to 
fall naturally with the c-axis oriented vertically, the 22° halo is 
observed in the cloud chamber with the Rassvet device. 


5. Discussion and summary 

Measurements in the 1 2 May SUCCESS contrail show that the 
ice particles observed in the core were small, mostly 1 20 pm, 

and while they appeared spheroidal in shape, instrument resolu- 
tion prevented determination of plates and similar shapes that 
may appear to be nearly spherical. On the contrail periphery, 
where the relative humidity exceeded ice saturation by 20-40%, 
ice particles as large as 300-500 pm were observed Particle 
shapes in the periphery included unresolved spherical, irregular, 
columnar and bullet rosette. Ice crystals larger than about 200 
pm were typically bullet rosettes. It is postulated (Heymsfield et 
al. 1998) that the larger crystals fall out of the contrail periphery 
and become precipitation streamers. The shapes of ice particles 
in the contrail periphery were typical of those sampled in patches 
of cirrus observed in the area; however, the cirrus contained 
fewer small particles and more larger ice crystals, extending to 
sizes of 600 pm and occasionally 1 mm. 

The shapes and sizes of ice particles generated in the CSU 
cloud chamber were remarkably similar to those observed in the 
contrail and ambient cirrus. Measurements of phase functions in 
the cloud chamber suggest that for small randomly-oriented (< 50 
pm) ice particles, as the crystal departs more from a spherical 
shape, relatively more light is scattered in the 40° -> 160° region 
Volkovitskiy et al. (1980) also found that more light is scattered 
from ice particles than cloud drops in this range of angles, and 
they suggest that this is due to relatively weaker scattering for ice 
particles in the 2 -> 40° range. Our measurements suggest that 
light scattered in the 40° — > 160° range from randomly-oriented 
ice crystals > 50 pm has a weaker dependence on crystal type and 
that crystal size is the dominant factor. Of course, these conclu- 
sions are based on a limited data set, about 12,000 individual ice 
crystals generated in a cloud chamber, and more data are needed 
for additional verification. 


14258 and NASA FIRE NAS1-96015. Grants to NCAR include 
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Air Force Office of Scientific Research grant F49620-96-C-0024. 
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1. Introduction 

Cloud-radiative processes in the Arctic have a strong impact on the stability of 
the Arctic Ocean ice pack (Curry et al. 1993) and also have ramifications on the 
global energy budget (e.g., IPCC 1990). Clouds in the boundary layer are 
persistent during May through -September and strongly influence the melting rate 
of the pack ice (Curry et al. 1996). The cloudy boundary layer in the Arctic is low, 
optically thin, and increases the summertime melt rate of sea ice, since the 
lonqwave exceeds the shortwave cloud-radiative forcing at the surface^ a 
positive feedback scenario occurs as melt ponds and leads form and the surface 
albedo decreases. 

Prior to the 1997-1998 Surface HEat Budget of the Arctic (SHEBA) project 
(Perovich et al. 1999), and the Eirst 1SCCP (International Satellite Cloud 
Climatology Project) Regional Experiment Arctic Clouds Expenment (FIRE ACE) 
(Curry et al. 2000), there were relatively few aircraft studies of Arctic clouds. 
While many types of Arctic stratus clouds were investigated during SHEBA, 
including clouds with multiple layers that extended to above 500 mb and cirrus, 
thin stratus clouds within 1 km of the surface occurred most frequent y and were 
the clouds most often investigated by the in situ research aircraft. We focusojjr 
analysis in this paper on these clouds, and refer to them as cloudy layers (after 
the terminology psed by Gurry et al. 1988; 1996), or “boundary layer clouds 
(Jhich was the^ terminology used by scientists at the FIRE. ACE field project and 
used in Curry et al. 2000). While Arctic boundary layer clouds may not always be 
thermodynamically or mechanically linked to the underlying surface, as is 
conventionally assumed in the classic definition, the term has been adopted here 
to represent low stratus cloud layers in the Arctic. 

In addition to investigations of boundary layer clouds, microphysical 
measurements from a deep stratus cloud that extended from 2 Km to B Km are 
presented as well as data from a cirrus cloud. These data are included in order 
to provide a better representation of the variability that exists in the many types of 
clouds observed in the Arctic. Also, the microphysical properties revealed in 
these investigations could be used to validate remote measurements and 
process models. 
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Curry et al. (1996) summarize findings from previous investigations of 
summertime Arctic cloudy boundary layers. The salient features from Their study 
that relate to this research are summarized here. 

• Three types of summertime cloud boundary layers are identified: i) a 
stable boundary layer with thin, patchy stable clouds that may be found in 
multiple layers, ii) a stable boundary layer, often with fog at the 

that is topped by a cloud-topped mixed layer, and iii) a cloud-topped, well- 
mixed boundary layerThat extends to the surface. 

. Warm moist air flows from continental regions over the pack ice and 
condensation is induced initially by radiative and diffusional cooling to the 
colder surface and longwave radiation to space. 

• Clouds in the boundary layer are often well-mixed, with the mixing 
assumed to proceed downward from cloud top due to radiational C0( ^ ,n 9 
and overturning. The well-mixed layer may extend to the surface, in which 
case, surface effects maybave also contributed to cloud formation. 


• The microphysical properties of clouds associated with the Arctic 
boundary layer are varied and definitive trends are difficult to establish. 
This may be due to the relatively small microphysical data set due to a 
lack of measurements. Significant jce concentrations are generally 
observed at temperatures colder than -15° to -20° C (Jayaweera and 
Ohtake 1973; Curry et al. 1990). However, exceptions to this generality 
have been reported. Some ice has been observed in clouds at 
temperatures as warm as -4° C. A predominantly water cloud was 
observed (in wintertime) at a temperature of -32° C (Witte 1968) and 
Curry et al. (1997) report an all-ice cloud at -14° C. Mixed-phase clouds 
are often observed and the types of ice particles have not been weH 
documented, except Tor surface-based observations. Curry (1986) found 
a significant amount of drizzle associated with a large dispersion in the 

droplet spectra. 

The microphysical properties of FIRE.ACE clouds are j rtvest 1 galed here using 
data collected by the National Center Tor Atmospheric Research j P9CAR) C-130 
research aircraft during May and July 1998. The study includes for the first bme 
in the Arctic, a new particle imaging probe that provides high-definit on i digital 
images of cloud particles. Data from the cloud particle ir " a 9 e f CPI) are 
analyzed to separate water drops from ice particles and identify crystal habits, 
compute water and ice particle size distributions. The CPI is i descntrec I in more 
detail, along with other microphysical instrumentation used in this study, in 

Section 3. 

The data presented in this paper are organized in the following way. A table 
gives an overview of the major physical features of clouds observed during the 
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16 C-130 missions in FIRE. ACE, focusing on clouds associated with the 
boundary layer. An example of a boundary layer cloud that is m °st!y adiabatic 
and homogeneous with a monomodal drop size distribution is discusse ^ 
Measurements of cloud liquid water content (LWC) are compared with th 
theoretical adiabatic values in the boundary layer clouds that are identified 
“adabatlf An example oT a Poundary layer 

actively mixing at cloud top and inhomogeneous wrth bimodal ^ 
distribution is presented. Time series measurements of cloud LWC ar °P'® 1 
concentration temperature and vertical velocity at different levels in six boundary 
laver clouds are discussed. An example of a boundary layer cloud with highly 
va^Sble hydrometeor fields is discussed in some detail. Images of P artic |®® ^ 
ice/water particle size distributions in a deep stratus cloud extending from i 2 km to 
6 tr are presented. Lastly, an example of the inhomogeneous (clump/) 
distribution of particles in a cirrus cloud is discussed. 

2 ‘ C urry ^et* a I ? '(2*000 ) °descri be the SHEBA/FIRE.ACE project in detail. Some 
salient features of the project that are pertinent to the C-130 flights are excepted 
here The main goal of the experiment was to examine the effects of clouds on 
radfation exchange between the surface, atmosphere and space, and to study 
how he surface influences the evolution of the cloudy boundary layer Data 
SclS during the field phase of the project are being used to evaluate and 
improve climate model parameterizations of Arctic cloud and radiation processes, 
satellite remote sensing of cloud and surface characteristics, and understanding 
of cloud-radiation feedbacks in the Arctic. 

The location and timing of the FIRE Arctic Clouds Experiment were 
determinedly the scheduled operations of The SHEBA W" £££ 
Beaufort Sea during October 1997 through October J?®?: . ,^ e on 

Ccastguard icebreater^^s GraMiWers was ^Ployed mult, -year £jo. ^ 

say? r jsucsr jst s 

during May, and nearly three hours in June, leaving 2 - 4 hours of on-station 
collection. 


3. Instrumentation 

ThP caoabilities of the NCAR C-130 and instrumentation on the research 
aircraft Curry et a. (2000). Of particular interest to this study 
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are microphysical instruments used to measure cloud particle charactenstics and 
cloud liquid watsr contont (LWC), including 

. Two King hot-wire LWC devices (King et al. 1978) manufactured by PMS 
with modifications to the electronics by NCAR, 

• a Gerber Scientific particulate volume monitor (PVM-100A) described by 
Gerber et al. (1994), 

. a Particle Measuring Systems (PMS) Forward Scattering Spectrometer 

Probe (FSSP-1 00), described by Knollenberg (1981), 

• a PMS 260X one-dimensional optical array probe with 64 photo diodes a 

10 [im pixel resolution, (Knollenberg 1981), 

• very limited use of the shadow-or concentration from the PMS 2D-C two- 

dimensional optical array probe (Knollenberg 1981), and 

• a CPI described briefly by Lawson et al. (1998), Korolev et al. (1999) and 
in more detail by Lawson (1997) and Lawson and Jensen (1998). 

The King probe is a hot-wire device with airflow characteristics that are 
theoretically predictable. The response of the probe has been shown to roll-off 
for drops larger than about 50 urn (Biter et al. 1987). The King probes are 
mounted nea 9 r each wing tip and close to the leading edge of the wing itself. 
Studies of the airflow around the probe installation on the C-130 have not been 
conducted Laursen (1998) points out that the location close to the wing leading 
e dqe was a matter for concern, and one King probe was re-located and extended 
from the wing for one mission. However, she reports that *h|s a P^are^oma)<e 
no discernable difference in the measurements. The King dry-ar term i is a 
function of the (local) true airspeed, and the dry-air term is subtracted from the in- 
doud°reading to obtain LWC It is therefore important to know the corred true 
airspeed at the probe location, or to subtract the dry air term measured outs'de of 
cloud at the same airspeed and air density. In this study, the King pro e 
were used as supplied by NCAR; the dry-air term is computed assuming the 
airspeed at the probe location is the same as that measured by the aircraft data 
system. Magnitudes of the errors that may be induced as a result of neglecting 
airflow considerations are undetermined. 

The PVM-100A is an optical device that measures LWC, drop surface area 
and effective radius at a sample rate of 1 KHz. The light scattered in the forward 
direction by an ensemble of drops is optically weighted and s ^ mn ]® d ° n ^ 
photodetector. The measurements presented here were averaged to 1 Hzl by the 
Mr AR data Drocessing routine. The probe was initially supplied with an 

was modified after the field expenmen. by Gerber 

Scientific Incorporated. The re-calibrated data have been used in this paper. 

Icing tunnel tests were conducted after the field experiment to 
LWC devices. The King probes performed well in the tunnel after the 9 a ‘^ ° 
of the probes was adjusted in software to compensate for an improper hardware 



5 


setting When the drop median volume diameter (MVD) was <~25 pm both 
King probes generally measured LWC to within 10% of the tunnel value (Strapp 
personal communication). The Gerber PVM probe on the » 0130 was no 
functioning for the icing tunnel tests. The University of Washington (UW) 
operated an identical PVM probe during FIRE. ACE and its performance in the 
icing tunnel tests is discussed by Garrett and Hobbs (1999). They report that 
during the small droplet tests the PVM underestimated LWC by 12/o in the 0 to 
0 75 g nrf 3 range. An analysis of the response of the UW PVM-100A to large 
droplet sizes showed that it had diminished response to droplets with diameter 
larger than 20 pm. At an TWO oTT2 pm the PVM underestimated spray LWC by 
-40% (compared to the 1 8% underestimation at a MVD of 30 pm reported by 

Gerber et al. (1994)). 

The FSSP has been under considerable scrutiny since its introduction into the 
field. Processing of FSSP data byHCARforthis project included partial recovery 
of losses due to coincidences and probe dead time (Baumgardner et al. 985), 
and adjustment of channel widths to account for airspeed corrections to the 
electronics (Cerni 1983; Dye and Baumgardner 1984; Baumgardner 1987^ 
Corrections have not been made for problems in the droplet spectrum due 
coincidences, as discussed by Cooper (1988) and Brenguier (1989), or due to 
laser beam inhomogeneities (Baumgardner and Spowart 1990; Wend| sch 1 998)^ 
The principle measurement of the FSSP is drop size and signals were sorted into 
15 equal drop-size bins in the 3-45 pm range. The measurements were summed 
and recorded every 0.1 s. The dynamic accuracy of the FSSP in measuring LWC 
and drop size is difficult to quantify. Baumgardner (1983) suggests that FSSP 
measurements of drop size are accurate to 17% and LWC is accurate to wthin 
34% However subsequent evaluations of the FSSP have shown additional 

potential error terms due to coincidences (e.g Cooper J988). 
the laser beam and effects of airspeed (e.g„ Wendisch et al. 1996). In addition, 
the accuracy of the FSSP appears to depend on factors that are not always 
quantifiable such as field calibrations, optical contamination, airflow effects due 
to position on the aircraft, etc. 


The PMS 260X probe sizes particles in 64 channels from 10 pm to 640 pm. 
According to Laursen (1998), 260X data processed by NCAR were filtered to 
eliminate the first three channels, due to changes in alignment caused by 
temperature variations. Thus, the first usable size channel started at 40 prrr 
Also, the end element photo diode voltages were monitored for indications that 
the probe was out of optical alignment and data were erroneous. 


Laursen (1998) also reported that concentrations from the (25 pm) 2D-C and 
(200 pm) 2D-P image data were suspect and could be in error by a factor of two_ 
Consequently, only the 2D-C and 2D-P shadow-or counts were output to the data 
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archives. The shadow-or particle concentrations reported in this paper were 
obtained by dividing the total counts by the true airspeed. Because there is no 
image information in the shadow-or measurements, conventional image 
processing routines that remove “artifacts” (i.e., water streaking on the mirrors 
and particles that break on the probe tips) can not be employed and sample 
volume corrections for small particles cannot be applied. Additional information 
describing the C-130 data contained in the FIRE.ACE can be found in Laursen 

(1998). 


The CPI is a relatively new instrument that utilizes innovative technology to 
record high-definition digital images of cloud particles, and measure particle size, 
shape and concentration. The high quality of CPI images supports the 
generation of individual size distributions for (spherical) water drops and (non- 
spherical) ice particles. The CPI casts an image of a particle on a solid-state 
one million pixel CCD camera by freezing the motion of the particle using a 25 ns 
pulsed high-power laser diode. As shown in the conceptual drawing in Figure 1 , 
a particle detection system (PDS) with upstream lasers precisely defines the 
focal plane so that at least one particle in the image is almost always in focus . 
Processing of CPI images significantly reduces out-of-focus siang ©rr^ fliat 
have plagued the conventional 2D imaging probes (Korolev et al. 1998). Each 
pixel in the array is 2.3 pm, so particles from about 10 pm to 2mm are imaged. A 
video-processing engine identifies and sizes particles within the one million pixe 
array saving only the regions around each image. The CCD camera can run at 
rates’ up to 40 Hz, and since more than a hundred particles per frame can be 
processed, data rates in excess of 1000 imaged particles s are obtainable. 

CPI particle concentration can be computed using a number of methods. 
Here we describe two nearly independent methodologies and an additional 
method that uses the 150 - 500 pm region of the PMS 260X particle size 
distribution to scale the CPI measurements. 

When the CCD camera is being read, the PDS continues to count parties, 
providing a continuous measurement of total particle concentration (called the 
'total strobes’ measurement of concentration). The total strobes concentrat'on 
a function of the threshold level of the PDS, which is user selectab e. A check on 
the effectiveness of the PDS can be obtained by momtonng the number of 
(‘valid’) camera frames that actually contain particles. Typify, the PDS 
threshold is set so that the percentage of valid frames is between 60 100 A. 

Since larger particles will more reliably trigger the PDS, there is a roll-off in the 
particle detection efficiency that starts at about 25 pm and depends on the PDS 
threshold setting. Thus, the small end of particle distnbutions that are narrow, 
such as a typical distribution of cloud drops, will be undercounted. 

In addition to the in-focus particle that triggers the PDS, high concentrations 
of small particles in the viewing area may produce additional images that are in 
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various degree of focus. In this case, the viewing area contains a snapshot of a 
known volume of particles, and an independent calculation of particle 
concentration, called the local concentration in this manuscript, can be 
determined by counting the number of imaged particles, and dividing by the 
viewing volume. Hoover, the optical depth of field (DOF), and consequent 
the image size and sample volume of the CPI, are a function of particle size (see 
Lawson and Cormack 1995, Korolev et al. 1998). The CPI particle ’ 
distributions shown in this paper have been made using size and sam ple volume 
corrections that are done on a particle-by-particle basis. Th^e correct^ns to 
the particle size distribution have a major effect on the concentration of small ( 
-100 pm) particles. This is possible because CPI images which c °"tain 256 
grey levels, can be used to accurately determine the fo< ^ of the image and th 
distance from the object plane (see Lawson and Cormack 1995). 

Due to varying DOF, the imaging sample volume of the CPI varies from about 
from 0 002 to 0.2 cm 3 . Local particle concentrations are a true measure wi 

the small volume that is imaged, but they need to be averag ® d 0 ^ r ^ics This 
or hundreds of snapshots to avoid large errors due to sampling statu ^ 
local concentration is still biased towards the small volumes v 5? 
concentration is highest, because the likelihood of tnggenng is higher where 
there are more particles. In regions with low particle concentrations, or ig y 
inhomogeneous cloud, the local concentration can be many orders of magnitude 
areater than the average concentration. Thus, the CPI size distributions mus 
Lied to the average concentration. In this study the 2M)Xa»d^ibution 
where it overlaps the CPI size distribution, was used as the first choice > for 
scaling (since the PMS 2D-C particle size distributions were not . avs aHable from 

the data archives). Where the 260X data were unre ^f 'f;J he / ^ " 
distribution was used. This method attempts to make use of the best 
measurement characteristics of each instrument. As a i sti • ’ 
the CPI total strobes average concentration may beuse ^° r ^f^ 1 00 q ,-h and 
most applicable when particle concentrations are relatively low (< = -1000 L . ) and 
uniform, and the CPI PDS threshold setfings are relatively low and unchanging. 

Due to the high resolution of CPI images and the 256 gray levels, it is 
possible to distinguish spherical from non-spherical particles, depending on the 
leveTof focus and size of .he particle. Generally, partides M are ngood focus 
and > 50 pm can be distinguished as spherical or non-sph«ncal This 
for separating water drops and ice particles in 

Darticles will generally grow to recognizable non-spherical shapes in less tha 
min" mixed-phase cloud. In this study we used 
automatically reject images that were not in sharp focus and then they were 
classified by another software algorithm that measures the roundless of he 
image We also classified several hundreds of particles by eye o v ®^ 
accuracy of the automated algorithm. The agreement between the automated 
and manual techniques was very good for images >~40 urn. In re 9 ,on ®^®'® 
classification of images < ~40 nm was essential, such as in regions wi 
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multiplication, a manual classification of particles was implemented. If particles 
smaller than about 20 pm could not be confidently classified, they were classified 
as water drops. 


4. Arctic Boundary Layer Clouds 


4.1 Overview 

Sixteen research missions were conducted during FIRE. ACE by the C-130. 
Half of the 16 missions were flown in May 1998 and the remainder were flown in 
July 1998 Of the 16 missions, we identified 1 1 cases when there were boundary 
layer clouds. Here we classify a boundary layer cloud as one that is well-mixed 
from the surface to cloud base, or if it has a cloud base that is < 3°0 m T able1 
shows the flight number, date, temperature, cloud base cloud top cloud 
classification, mixing depth and some general microphys.ral^^^^^^^ 
each of the 1 6 missions. Upper air soundings were taken from the SHEBA sh p 
approximately every six hours. The mixing depth was determined by analysis of 
U^^imding^ that was closest in time to the flight profile, so some of the fine 
atmospheric structure may be missed in this analysis. Even so, here appears to 
be a definitive difference in the subcloud properties in May and July. In May 
out of eight cases had boundary layer clouds and all six clouds were mixed from 
the surface to ctoud base. The depths of the mixing layers in May ranged from 
150 to 1200 m. In July, five of the eight cases had boundary layer clouds, but 
none of these were mixed from the surface to cloud base. 

4.2 Mostly Adiabatic, Homogeneous Water Clouds 

The C-130 flew a total of 21 slant profiles on 15. 18, 27 May and 29 i July 
through boundary layer clouds that consisted almost entirely of liquid water n 
cloud 9 physics literature, adiabatic generally implies both adiabatic and I dosed in 
classical thermodynamics terms. That is, when a parcel of a.r undergoes an 
adiabatic ascent from cloud base, it does not exchange mass °f h “* ™‘ h 
environment. Under these conditions the LWC and temperature ran ^e 
predicted theoretically from the cloud base temperature and pressure (see 
Lawson and Blyth 1998). In addition, the cloud drop SIZ ® °fg 7 gf 

estimated usinq the droplet growth equation (see Pruppacher and Klett 1978). 
523? during adiabatic ascent the mean drop size will grow monoton, cally 
and the^drop size distribution will remain monomodal. Bimodal droplet size 
distributions P (not observed in these homogeneous clouds) often result from 
^7* environmental air (above cloud top), followed by ascent and 
activation of new CCN (Jiang et al. 2000, present modeling results on the effects 
of CCN entrainment). In a well-mixed boundary layer the vertical profile of 
temperature and liquid water content are adiabatic regardless of whether parcels 
SXSTnJS. with their environment. This is because the parcel and 

environment are identical. 
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Figure 2 shows an example of temperature, LWC and drop size distribution 
measurements from a boundary layer cloud observed on 18 May 1998 that is 
mostly adiabatic and homogeneous. In the example shown in Figure 2, the 
temperature profile is nearly adiabatic, the profile of LWC measurements from 
the King probes and PVM are within 20% of adiabatic. (The FSSP LWC is about 
twice the adiabatic value and is believed to be in error, as discussed in the next 
section.) The FSSP drop size distribution is monomodal, the profile of mean 
droplet size increases monotonicaTly with height and the dispersion of the droplet 
spectrum is constant with height. There was no drizzle detected. All of these 
factors support the conclusion that this boundary layer cloud observed from 2208 
-2210 UTC by the C-130 on 18 May T998 was mostly adiabatic. 

CPI images provide another method of investigating cloud homogeneity. As 
seen in Figure 1, when a particle is detected by the PDS and the CPI imaging 
laser is fired, the imaged volume may contain more than one particle, particularly 
if the particle concentration is relatively high. In a homogeneous cloud, the effect 
of sampling random volumes of equal size would result in a Poisson distribution 
for the number of particles per sample. Figure 3a shows the Poisson distribution 
for the number of particles per sample with the same mean as that observed by 
the CPI from 221010 - 221030 UTC during the slant profile shown in Figure 2. 
Also shown in Figure 3a is the distribution of the number of particles per image 
frame actually measured by the CPI. The comparison in Figure 3a shows that the 
CPI measurements closely follow the theoretical Poisson distribution. 

A more sensitive test for inhomogeneity can be made by conditional 
examination of the particle size distributions for the cloud region shown in 
Figures 3a. One conditional distribution was made of the sizes of particles 
imaged in CPI image frames with three or less particles, while the other was 
made from the sizes of particles imaged in frames with five or more particles. If 
the cloud were homogeneous, the conditional size distributions should be 
identical The data in Figure 3b show that the conditional spectra are nearly 
identical, suggesting that this region of cloud was quite homogeneous. 

The measurements shown in Figures 2 and 3 are representative of 
FIRE ACE boundary layer clouds with adiabatic LWC profiles. However, even in 
these cases there is evidence that some entrainment and mixing is occurring in 
some regions of the cloud. In Figure 2, clear air holes are seen that extend to at 
least 100 m below cloud top. The CPI also confirms the lack of cloud in these 
spots. Presumably, these are entrainment events that bring clear air from above 
the cloud into the cloud layer. Curry (1986) also observed that cloud top air 
penetrated downward into the cloud to a depth of at least 50 m. 
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4.3 Evaluation of LWC Instrumentation in Adiabatic Clouds 

In addition to providing an excellent natural outdoor laboratory for studying 
mixing and evolution of the droplet spectra, these Arctic clouds provide a reliable 
test bed for assessing the performance of some of the microphysical 
instrumentation. The LWC profile and temperature are predictable in regions of 
cloud that have ascended adiabatically from cloud base to the observation level 
(e.q. Jensen et al. 1985; Lawson and Blyth 1998). The cloud base temperature 
and pressure were measured at the point where the FSSP droplet concentra ion 
exceeded zero on ascents, and when the concentration dropped to zero on 
descents. A zero FSSP concentration threshold was used because, unlike some 
environments where there are relatively high concentrations of large aerosols 
that produce occasional counts, the FSSP almost always read zero in clear air. 
This method of computing adiabatic LWC is not only the simplest to implement, 
but arguably the most accurate, because the subsequent analysis depends 
primarily on relative changes in temperature and pressure. Static pressure 
measurements on research aircraft are generally felt to be accurate to 1 mb. 
dry air an absolute uncertainty of about 0.3° C is expected in the Rosemoun 
temperature measurements (Lawson and Cooper 1990). Lawson and Cooper 
(1990) show that an additional error from sensor wetting of < ~0.5° C is expecte 
when the LWC is < ~0.5 g m' 3 , as was the case in these clouds. No evidence of 
temperature spikes, typically found when the Rosemount temperature probe 
element gets wet (Lawson and Cooper 1990), was observed on cloud exit For a 
1 ° C error in measurement of temperature and a 1 mb error in pressure at doud 
base there is only a 0.025 g m' 3 maximum error in adiabatic LWC for conditions 
typical of FI RE. ACE boundary layer clouds. Thus, instrumentation measurement 
errors in pressure and temperature at cloud base can be expected to have an 
insignificant contribution to the determination of adiabatic LWC. 

Figure 4 shows scatterplots of LWC measurements versus adiabatic LWC for 
all of the (nine) ascents and descents through clouds that were determined to be 
close to adiabatic, although as shown in Figure 2, 'holes’ were ocras'on^y 
observed in these clouds. Since there is no reason to expect su P er 3 d '^^ 
LWC in these clouds, and since there is evidence of some mixing almost 
everywhere the calculated adiabatic value is the maximum that the 
measurements should reach (except for the contribution from measur ®"®" 
errors). The data Hi the scatterplots show that both King probes generaUy did not 
exceed the adiabatic value and were within 75% of adiabatic LWC (except in 
regions with obvious ‘holes’). The PVM scattered around the adiabatic value and 
sometimes exceeded it by up to about 35%. The FSSP generally exceeded the 
adiabatic LWC value by up to a factor of two. 

While instrumentation errors were shown to have a negligibly small effect in 
determination of adiabatic LWC, fluctuations in the actual cloud base level could 
potentially introduce significant errors in adiabatic LWC. For example, during 
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slant ascent or descent through a stratus cloud, the aircraft may enter cloud base 
at one point, and when it reaches cloud top, the air parcel it samples may have 
originated at a different cloud base elevation. Significant errors in adiabatic LWC 
are more likely when there are noticeable undulations in cloud base elevatorr 
Figure 5 shows LWC measurements from one of the nine adiabatic profiles and 
the effect on LWC for different cloud base elevations. It can be seen that 
deviation from the measured value of 1001 mb to about 1013 mb (a clrfferonco of 
about 120 m) would be required to produce an adiabatic LWC that agrees with 
fhe peaks in the FSSP measurements. It is veiy unlikely that variations of this 
magnitude existed in the actual cloud base elevations in this case, ^ evidenced 
by the measurements of FSSP and PVM LWC also shown in Figure 5 when the 
C-130 exited cloud and flew level just underneath cloud base for about 15 s. If a 
fower doud base existed, it would likely have been detected during th.s honzo^ 
,eo in addition after exiting cloud, the C-130 flew for about two min just below 
the cloud base elevation before climbing for another cloud penetration. During 
this time there was no indication of lower cloud base, and the cloud b 
measured upon ascent was within two mb of the previous measurement The 
nine slant profiles used to generate the data in Figure 4 where carefully select , 
so that errors in determination of cloud base were minimized. 

Thus we can conclude that for C-130 measurements from this field project^) 
the King probes rarely exceeded the adiabatic value and were usually 
to within 20%.'. 2) The PVM measurements scattered around the adiabab 
value with an occasional apparent tendency to overestimate the adiabatic • 
3) The FSSP measurements are heavily biased toward values that exce eded t 
adiabatic LWC It is highly unlikely that the FSSP measurements could be 
“ urtderestimates 9 in y .he adiabatic LWC, since .hist mid mean the 
aircraft consistently measured a cloud base that was about 10 mb too high, 
that the King and PVM probes consistently underestimated LWC. 

The FSSP counts and sizes droplets. Errors in ^zin^re rased to the third 
power and can result in significant errors in LWC. The NCAR 
calibrated with glass beads in the field, but even so, stzing i errors can occur 
fBaumaardner et al 1990). Errors from undercounting in high (> 600 cm ) drop 
SSns are due to coincidence, which also artificially 
size spectra (Cooper 1988, Brenguier 1989). Coincidence artificially ' ncrea ®®s 
LWC because two small drops are erroneously recorded as one larger drop a 
since LWC is proportional to the third moment, the undercounting IS 
by the increase in drop size. Hoover, (For 


cm' 


e increase in arop size, nuwevc, , y"..“ V WFRA <*hio fFor 

except perhaps downwind of the effluent plume from the SHEBA ship (for 


1 Shortlv after this manuscript was accepted, measurements of the dimensions of the ^ng probe 
senOT wires used talhe FIRE.ACE ptojeo. showed that the area of the rXe 
about 20%. in the calculation of LWC, sensor area if nnv«g » 'Pg^^-^^ignTacility 
FIRE.ACE measurements would increase by about 20%. The NCAR Re^arch Av.at.onr 
/RAF) plans to publish a correction that can be applied to the archived FIRE.ACE data (K. 
Laursen 2000 - Personal Communication). The end result will be that the King values 

measurements will be increased by about 20% and be in closer agreement w.th the PVM values. 
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i j, ,riri/i cupra epp Yum and Hudson 2000 and Pinto ©t 
"T™ 3^1^* problematical. Thus, there is 

^rrlTyno e^"for the apparent overestimate of LWC observed ,n the 
FSSP measurements. 

Reoardless of the apparent error in FSSP measurement of LWC, the primary 
JS lhe FSSP isTo measure drop size distnbuboa There - co - 
determination of how much of the apparent am mJTatow 

E5&5E “^efU and are 

included frequently in this paper. 

The FSSP displayed the most stable baseline outside of cloud and was the 

» »-»«■» r. rr? 

on every flight. 

A qualitative investigation of the effects of ice on measurements of ^LWC j was 
undertaken by comparing probe responses in a cate d°u I ™ h^l 'ce part,c ^ 
Figure 6 shows the responses of the King, PVM and FSSP ^obes to an ^ 
cinus cloud along with ice particle concentration and IWC measurements 
derived from CPI measurements. There are undefined uncertainties '" ’ he 
fwr measurements however, using the relative responses of the LWC 

S £ SK ZteZSZ M mte 

by the small particles. 

4.3 Non-Adiabaltic, Inhomogeneous Clouds 

As previously discussed, nine of the 21 slant profiles flown through the .nearly 

and actively mixing. 
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Fiqure 7 shows a C-1 30 slant profile through a boundary layer cloud on 29 
July 1998 Both LWC and temperature display a systematic trend to deviate 
increasing more from adiabatic values with increasing distance above cloud 
base, which strongly suggests an active mixing process at cloud top (Zhang etaL 
2000 present modeling results on cloud top mixing). At cloud top, 
temperature fluctuates strongly and is often warmer than adiabatic, due o 
entrainment of relatively warm air from the temperature inversion above cloud 
top The droplet spectra are bimodal at cloud top and drizzle was encountered 
lower in the cloud. This is in contrast to boundary layer clouds where the LWC 
and temperature profiles were adiabatic and the drop spectra were monomodal 
(i.e. Figure 2), where drizzle was not observed anywhere along the profile. 

When the LWC and temperature profiles were significantly sub-adiabatic in 
FIRE ACE boundary layer clouds, the droplet spectra near cloud top were o en 
bimodal and drizzle was observed lower in the cloud. It was not a necessary 
"on for the drop spectra to be bimodal to observe drizzle, since in some 
subadiabatic clouds that were actively mixing, drizzle was occasionally °bs®rv®c! 
at cloud base without concurrent measurements of bimodal drop spectra in th 
cloud However, it should be remembered that the aircraft samples a relative y 
small volume of cloud, and bimodal drop spectra could be present at locations in 
the cloud not sampled by the aircraft. 

Figure 8 shows another profile flown about 7 minutes earlier than the slant 
profile in the same boundary layer cloud shown in Figure 7. Basically, the same 
large-scale observation seen in Figure 7, which is that the boundary layer doud 
is actively mixing at cloud top, is seen in Figure 8. The profile shovm in Figure _8 
however has regions where the C-1 30 flew several-minutes-long constant 
altitude legs providing an opportunity to examine cloud inhomogeneity using the 
conditional spectra technique introduced in Section 4.2 and shown for a relatively 
homogeneous cloud region in Figure 3. The presentations of theoretical Poisson 
and CPI measured number of particles per image frame and the conditional drop 
size distributions, analogous to those shown in Figure 3, are shown above the 
time series measurements in Figure 8. The measured number of particles per 
image frame follow the Poisson distribution and the co ^ l, ' onal 7 d Q r0 P ““ 
distributions are nearly identical for the doud regions at 180 and 270 miTist 
This implies that these regions are well-mixed and there is not a toga 'degree _of 
inhomogeneity in the lower and middle levels of this boundary layer doud. On 
the otter hand near cloud top at 450 m msl, a larger degree on inhomogeneity 
be seen in the temperature and LWC measurements, andthis is refleded in 
the conditional drop size distribution measurements, which are noticearay 
different. It is interesting to note that the trend in the ^ditio^l spwtra is for the 
smaller droplets to be found in locally higher concentrations than the larger 
droplets. Just the opposite mighf be expected in an actively mixing cloud. 

It is also interesting to note that the drop size distribution is distinctly bimodal 
in the mid-level region, yet the conditional spectra revealed no inhomogeneity. 
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This suggests a process whereby the mixtures formed at cloud top (coo y 
evaporation) descended toward mid-cloud and continued mixing while 
descending, which resulted in a relatively well mixed region. 

Perhaps the most striking aspect of boundary layer d™* *sehred i during 
FIRE ACE as revealed by this analysis, is how a group of clouds that are fairly 

similar in visual appearance can exhibit such variabH.ty .n micro^hys^ 
Drooerties In order to investigate further aspects of the FIRE.ACE boundary 
layer clouds, we generated time series measurements of key microphysical, 
thermodynamic and dynamic parameters for six cases. 

Fiaures 9 and 10 show time series measurements of LWC, temperature, 
vertical velocity FSSP concentration, altitude and the standard deviation for all o 
these quantities, except for altitude. The data in these figures are organized to 
show examples of the microphysical variability between these boundary layer 
%Zs ^ There ^bns selected for display in Figures 9 and 10 ^rechosenwhen 
the C-130 flew level for several seconds, and at different altl * ud ® s in i the 
boundary layer clouds (the only exception being 15 May in Figure 9 when there 
were no^representative level regions). The purpose of selecting these regions 
was to see if there were any consistently recognizable features that varied as a 
m of altitude in the boundary layer clouds. Figure 9 shows f e^mptes of 
cloud regions (also shown in Figures 2 and 4) that were m°s«y ad '^ " ,h 
relatively constant FSSP drop concentration and temperature measi uremei nt _ 
These clouds displayed monomodal droplet spectra and U l® r ®. n m 
observations of drizzle at the time of the in situ measuremerts Figi^^^ 
other hand, shows cloud regions where mixing produced sub -® d '® d ®'' ® 
variations in the FSSP drop concentration and temperature, and often bimo 
droplet spectra and drizzle. 


The data in Figure 9 show that there is very little change with altitude in the 
thermodynamic (i.e., temperature), dynamical (i.e„ vertical velocity) and 
Sphysical parameters (i.e., LWC and FSSP concentration) in these aH-v^ter 
clouds As shown in Figure 3, the droplet spectra were monomodal and no 
drizzle' was observed by the C-130. As previously stated, this is Typical Id f the 
“mostly adiabatic” boundary layer clouds. The values of standard deviabon (a) 
in Figure 9 are mostly constant, except for the bumps in drop concentration a 
LWC^ (resulting from the 20 s over which a is computed) around the region where 
the ‘holes’ are observed in the 1 8TS7lay cloud. 

Figure 10 shows four examples of profiles in boundary layer clouds i that were 
actively mixing. The data in this figure show that fuduations _in ‘®^ratufe_ 
LWC and FSSP drop concentration tend to increase as the C-1 9 . 

cloud top This is consistent with the hypothesis put forth previously that these 
boundarv layer clouds are mixing from the cloud top downward. The 
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measurements in Figure 10 show that there is not an increase in the standard 
deviation of vertical velocity with height in these boundary layer clouds. The air 
motion system in the C-130 contains a residual offset that is on the order of 0.5 
m s' 1 , so the small offset in vertical velocity in Figures 9 and 10 should be 
ignored However the small-scale fluctuations are considered to be accurate to 
within a few tenths m s*\ so the consistent lack of increase in the fluctuation of 
vertical velocity with height is significant. The boundary layer clouds depicted 
by the measurements in Figure 10 all had bimodal drop distributions at cloud top 
and drizzle was observed below cloud. 

Generally, in Figures 9 and 10 there was a lack of correlation in vertical 
velocity with FSSP drop concentration, and there is a correlation between LWC 
and drop concentration. An exception to the above generalizations is seen in 
Fiqure 9 on 18 May, where the drop concentration increases from the nominal 
level of about 80 cm -3 to 140 cm^ in three regions, and the LWC remains fairly 
steady. In these regions the vertical velocity shows a significant increase. The 
FSSP mean drop size (not shown) was anti-correlated with drop concentration, 
which accounts for the lack of increase in LWC, and also agrees with the 
observation that a localized updraft could have activated more CCN. 

The test for inhomogeneity using CPI images, described in Section 4.2 and 
shown in Figures 3 and 8, was applied to the region in Figure 9 on 18 May where 
there is a subtle, but noticeable variation in drop concentration. Figure 1 1 shows 
the results in the same format as Figure 3. which is from the same cloud profile 
as shown in Figure 1 1 , but about 50 s later when there was no obvious variation 
in droplet concentration. The conditional spectra in Figure 11 are slightly 
separate suggesting that the region has a detectable inhomogeneity. e 
results of the inhomogeneity test using CPI image data shown in Figures 3 and 
1 1 demonstrate the sensitivity of this type of test, but do not provide information 
on the scale of the inhomogeneity. However, in this case, the large scale 
structure that can be seen in the time series measurements of Figure 9, likely 
caused the separation seen in Figure 11. The tendency for smaller droplets to 
be observed in higher concentration regions is consistent with the time-senes 

measurements. 


4.4 Mixed Phase Cloud 

Fiaure 1 2 shows an example of time-series measurements from a boundary 
layer cloud that is fairly unique in the C-130 FI RE. ACE data set. This was a 
case where cloud base was higher (640 m MSL) than the clouds previously 
discussed and the entire sub-cloud layer was well-mixed to the surface (see 
Fiqure 12). This was the first C-130 mission (4 May), the atmosphere was still 
transitioning from the winter regime, and with the higher cloud base this 
boundary layer cloud was colder than the other cases discussed. The trends in 
the measurements in Figure 12 differ from Figures 9 and 10 in that the standard 
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deviation in vertical velocity is noticeably larger ,ower J^ th f x ^J° ud ’ and r ^‘ S 
less near cloud top than near cloud base. Even though the LWC measurements 
are significantly less than adiabatic everywhere in this cloud there is no 
of active mixing near cloud top and the drop spectrum (not shown) at cloud top is 

monomodal. 

The most striking feature of the boundary layer cloud observed on 4 May is 
the inhomogeneity in the hydrometeor T.elds. Figure 13 shows « 
flight track when the C-1 30 was descending and making passes over the SHEBA 
ship from cloud top (1025 m) down to cloud base (690 m) and examples of CPI 
images water drop and ice particle size distributions during this time period^ 
The 9 data in the figure show that the hydrometeor fields varied considerably over 
spatial distances®^ 10 Km horizontally and a few hundred meters vertcally^ 
When the C-1 30 skimmed cloud top and then turned and made a pass 3 
below top it encountered mostly supercooled (-25.5 to -22 C) cou r °P' e s 

with a mono-modal size distribution and only a few 100 - 500 pm ice 
The larger ice particles may have fallen from a higher cloud. The ratio of water 
drop to ice particle concentration is > 1000:1, so the region near cloud top was 
predominately composed of supercooled cloud droplets with a monomodal size 

distribution. 

On the next pass 75 m below cloud top, Figure 13 shows that the CPI water 
drop ^or^ntration had decreased by about a factor of six and the size 
distribution still had a mode of about 30 pm, but had also broadened sligh y. 
Some ice particles were observed, mostly unnmed with sizes to about 300 pm 
and an occasional rimed particle as large as 600 pm. However, just west of this 
pSsittenTd sS m lower (920 mb), while the C-1 30 was turning to reverse 
course, a section of drizzle about 10 km in horizontal dimension was observed 
with supercooled (-25° C) drops with sizes up to 180 pm. After the pocket of 
drizzle was transected, the C-1 30 made another pass over the ship while 
descending to 900 m and encountered supercooled drops with a mode of abou 
30 pm and ice particles that were mostly < 500 pm, with a few larger rimed 
narticles But then when the C-T30 reversed course and two minutes ater flew 
only 30 m lower ’ (870 m) over the same flight track, significantly higher 
concentrations of rimed ice particles from 400 - 800 pm and occasional graup 
particles up to 1 mm were observed. 

The degree of inhomogeneity in cloud microphysics could possiblybe due to 
mixing downward from cloud top or from the effects of seeding from higher 
clouds Even though significant concentrations of ice particles were not 
observed during the passes over cloud top, the potential for temporal and spai 
variability of this type of seeding precludes eliminating this poss.bihty , 
presence of the small pocket of drizzle at -25° C, however, cannot be direct y 
explained by either surface or elevated effects. Curry (1986) r ®P° rte r 
when there was a relatively large dispersion in the drop spectra In ^is case the 
FSSP drop spectra at cloud top were relatively narrow and mono-modal, in 



17 


contrast to the example shown in Figure 7, where the drop spectra were bimodal 

and drizzle was observed at -1° C. "However, the C-130 could 
investigate cloud top, so that bimodal drop spectra may have existed somewhere 
at cloud top. The cloud on 4 May existed for a long time and seeding from 
above was P sporadic. The curious pocket of drizzle that was produced by 
coalescence may have been an older region that was not seeded from above. 

5.0 Deep Stratus Cloud with Very High Ice Particle Concentrations 

Next we examine inhomogeneity in microphysical measurements collected in 
a deep Arctic stratus cloud. This is not a boundary layer .cloud according^ the 
definition used in this paper. The cloud extended from about 2 Km to 6 Km ms 
(+2 to -23 ° C) on 18 July 1998 and contained a myriad of 
characteristics. Figure 14 shows a vertical profile of K,n 9 
and 2D-C (shadow-or) particle concentration measurements along 
and ice particle size distributions. Representative examples of CPI images of the 

particles are also shown. 

The data from the microphysical measurements shown in F ‘9u re J4 willbe 
discussed from cloud top down to cloud base. In the upper portion of cloudy 
from cloud top at -23” C to -13 ” C, the FSSP droplet concentration varies from 
about 50 to 250 cm' 3 averaging about 125 cm" 3 , and the LWC ranges from zero 
To ifg m- 3 averaging about 0.1 g t* 3 . The FSSP drop size distributions 
shown in Figure 14 are relatively constant in this region, but not especially bro i, 
with 20 urn drops found in concentrations >1 cm- 3 only in the regions with LWC > 
~0 1 a m' 3 However, CPI images show a relatively low concentration of drizzle 

drops up to 125 pm in diameter near -19 ” C. The 260X 
of ~40 L' 1 and the 2D-C shadow-or registered a maximum of 1 L . This appears 
to be an example of “non-classical freezing drizzle” formation (Cober et al. 1996 
Lawson et al, 1998), where drizzle is formed through ^ 

drops usually near the tops of stratus layers. However, in the c ° b ® r ®* a '_ 
Lawson et al measurements, the FSSP drop spectra were generally broader and 
the 2D-C concentrations were significantly higher, in excess of 500 

The CPI images show that ice particles were very rare in the region from -23° 
C to -13 • C, and were mostly large (750 pm) heavily rimed partides. ^ The 
averaae 2D-C and 2D-P shadow-or particle concentrations are <01 L an 
support the observation that there was negligible ice in the layer. The lack of ice 
in the region from -23° C to -13 ° C is a curious aspect of this deep stratus chDud 
because^as will be shown, ice concentrations in the warmer regions below were 

much higher. 


C there is a noticeable drop in LWC to about <-0.011^ g 


rrf 


Near -12 ° v/ mwc io <=> nwuwww — r- •• .3 

reduction in FSSP total particle concentration to about 5 cm ana “ 
corresponding peak in 260X and 2D-C concentrations. Figure 14 shows CPI 
images of some examples of the crystals and Table 3 gives a breakdown of 
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visual classification of the particles. Figure 14 and Table 2 I show that * 

mixture of irregular crystals with abundant s.deplane i g owtl 
rrvstals/orauDel and small (< 50 pm) spheroidal 2 ice particles. The CPI P a ™ cle 

Figure 14 slTow that ice existed in 

2500 L 1 in the region around -12 ° C. Figure 15 shows the F SSP/CPI/260X 
combined particle size distribution in the ‘central region with g 
concentrations and at the ‘edges’ of this region (i.e., immediately above i and 

below as viewed in the C-130 slant ascent) where the ^ “2S?tO?cm- 
reduced to about 750 L' 1 and the FSSP concentration increases to about 100 cm 


The microphysical data in Figures 14, 15 and Table 3 indicates that the reg'on 
of hioh ice concentration is composed of ice particles whose population is 
dominated by small (< 50 pm) ice spheroids that are not round, but can also not 

be identified as vapor-grown crystals. They could be frozen 
were once frozen to larger ice particles, or ever .fragments oi ^larger ce parMee 
The laroer ice particles are mostly irregular in shape and often contain side plane 
JoJth InThe^acenr regions immediately outside of the ares > v«t thigh ice 

particle concentrations, the FSSP water drop «> nc ® ntra “°" 

inn rm‘ 3 comDared with 5 cm -3 ). The most striking feature in the aa J a ^® 

a: es droids 

compared to the ^central region. 

in the region from about -1 1 • C to - 7 • C the LWC flu^ates beWveen zero 
and 0.05 g rtf 3 and the FSSP droplet concentration ranges from 50 to 1 00 ^ 
CPI images show low concentrations of rimed ice pa ic P 
occasional drizzle drops in excess of 100 pm. The ar*vity -on .#■ ' «K and 2D C 
are also relatively low, indicating few particles larger than 40 to 50 pm. 

Very high concentrations (-4,000 L ') of ice particles ^re a lsoobsenjedjn 

Hallet (1974); Mossop 1985): 

i ) cloud temperatures between -2.5 and -8° C, 

ii ) droplets > 23 pm in concentrations^ 1 cm 

,ii ) relatively fast-falling (>0.2 to 5 m s 1 ) ice particles, 

The FSSP total particle concentration ranges from 2 to 20 cm 
shovre^articlesi^e dSnhubons^nd CPI images in a format similar to that shown 


2 Here spheroidal refers to image shapes that are not spherical (i.e., perfertlyroiind), but 
instead, have a slight but noticeable departure from a perfectly sphencal shape. 
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SSJJM 

central region. 

The Rosemount icing detector did not register any su P er “f te ^ 5 ' A f t ^ ^ 
regions with high ice particle concentrations a Mazin et al ' (2000) 

have^show^that^he^heoretira^'sens^^o^the^Rosemount icin ® 

Sfn^fot^ 

themodynarrHc con^tions'ht^-wiVe proves will register a false LWC signal that is 
about 15 to 20% of the equivalent ice mass content. 

ThP H-M conditions do appear to be satisfied in the regions adjacent to the 

thes^reg^ons 'are stHl"* abo^^hose^xpe^ 

about 25% of the ice concentration in the central region. 

The mechanism(s) that is producing the very high ice particle 
The ™ ecna ™ s ™^> o r a nd _ 4 ° c is not obvious. Rangno and Hobbs 
in the regions near —12 u ana cidc apf stratus cloud 

&r-£s 

particles in a format similar to Table « . 3 ^ n g to note that both analyzes 
classification is somewhat subjective, it s interesiing labeled frozen 

5i?as sssrsisfff aa 
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shattering of isolated drops during freezing and fragmentation of '^ c ^ lsby 
collisions They conclude that either mechanism may have been °P erat ' n 9- 
seemed to slightly favor the fragmentation process, since it appeared to better f 

the observations. 

Nearly all of the larger ice particles in the 'adjacent' ^ 

The individual drops frozen onto these large partides ( F igures 15 and 16) Ufa 
could be susceptible to breaking off. It is also interesting to note that the small 

spheroids in the 'central' regions look as if they could be Sith e V ohvsi^ 
broken from the large rimed crystals. There is no established physical 
mechanism to support this hypothesis, however, a similar process of rime 
breakup has been suggested by Vali (1980). 

There is a thin layer of supercooled cloud from 
the hiahest (0.22 g m -3 ) LWC and the highest (280 cm ) FSSP droplet 
concentration observed anywhere in the cloud. Also, the CPI detected no ice 
particles in this region and the 260X and 2D-C concentrations were near zero, 
suaaestinq that this was a cloud composed of small water drops. It is difficult 
exp^n why the large particles observed above in the cloud are not seen fal ing 
throuah this reqion of supercooled cloud. This is especially difficult to expla 
because drizzle was observed below the freezing level, presumably the result of 
rnemng of lS la^er ice particles. However, since the C-130 ,s climbing at a 
shallow angle sp that the data shown in Figure 14 do not actual y comprise a 
vertical profile, it is possible that extreme variability in horizontal structure in the 
cloud (such as seen in Figure 13) can explain this observation. 

6. Cirrus Cloud with Highly Inhomogeneous Particle Distribution 

Here we focus on the extreme spatial variability (“dumpiness”) seen in a 
cirrus cloud. Figure 17 shows a combined particle size d| s tr,bu ^" 

CPI and 260X data as the aircraft flew through cirrus cloud at 5400 m (-25 C) 
for aooroximately 25 minutes on 29 July 1998. The Rosemount icing detector 

fheparticles 'a^aswm^tobe^^'^'^e'cpi partSe 

scaled to the 260X data in the 150 to 500 pm size region, where the 260X 

rraasuremerts are felt to be most reliable. The three particle size distributions 

show relatively good agreement in '^ e re 9“ n ® ^p r ® C pf g^^gox are'shtwnln 

Tr 9 “ r S C?| d so y aled e toKK^e Oration 

measurements ha^ been scaled to the CP, total part* 
rnmouted from the composite size distribution in Figure 17. The time-series 
measurements show that there is considerable spatial structure to the cloud on 
scales greater than the 120 m resolution of the instruments. 

Tho ccqp io considered here to be a reliable measurement of the average 
concentration of srnalt partcles^n cirrus clouds. Previous reports in the literature 
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(eg., Gardiner and Hallett 1985) suggest that the FSSP is unreliable in the 
presence of ice. However, these measurements were in mixed-phase clouds 
where contributions from ice particles and water drops were difficult to separate. 
More recently, the literature contains FSSP measurements that are felt to be 
mostly reliable when the probe is measuring small ice particles in cirrus (e.g., 
Gayet et al. 1996; Poellot et al. 1999; Amott et al. 2000). The reasoning for 
using the FSSP data as a measure of average ice particle concentration cirrus 

follows: . . 

i.) CPI imagery shows that small (< 40 pm) particles in cirrus are largely 

spheroidal in shape, so even though the probe may not size these 
particles as accurately as it does water drops, it will still record a signal 
from the forward diffraction peak that will be registered as a count, 
jj ) The FSSP-100 has a 6 ps dead time, so particles traveling at the 
(-120 m s 1 ) airspeed of the C-130 in cirrus must be > -750 pm to 
produce a double count. In addition, Figure 17 shows that the number 
of particles larger Ihan 750 pm is four orders of magnitude less than 
the contribution of the smaller particles, so even if there were doub e 
counts from the large particles, they would have a negligible 
contribution. 


The time-series measurements from the FSSP, CPI and 260X agree very well 
in phase (i.e., the peaks and valleys line up well). The concentrations measured 
by the FSSP are about 30% greater than those measured by the CPI, and both 
the FSSP and CPI are more than an order of magnitude great :er than the 260X 
concentrations. This is largely because, as shown in Figure 17, both the FSSP 
and CPI data show large numbers of small particles, and high ^^® ntr ® t ‘° n 
small particles are more effectively counted by the FSSP and CPI. The first 
useable channel in the 260X started at 40 pm, compared with 3 pm for the 
FSSP, while the smallest resolvable CPI image is about 10 pm. 

We now return to the discussion of the highly inhomogeneous particle 
distribution, or “plumpiness", in this cirrus cloud, which is independent of the 
absolute magnitudes of the average particle concentrations shown in Figure lb. 
Figure 19a shows the distribution of the number of particles per CPI image frame 
and the Poisson distribution with the same mean. As explained in the text 
describing Figure 3, the Poisson distribution would be the result of random 
sampling^ a homogeneous cloud; the comparison in Figure 19a shows that the 

cloud is inhomogeneous. 

We now show that small particles are found in clumps with very high local 
concentrations that are interspersed with regions of ' arger 
concentrations. The conditional size distnbutions (analogous to F, 9ure 3b and 
11) are shown in 17b. One conditional distribution was made of the sizes of 
particles imaged in frames with only one particle, while the other was made rom 
the sizes of particles imaged in frames with seven or more particles. If the cloud 
were homogeneous, the conditional size distributions would be identical. 
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However, as seen in Figure 19b, the high concentration frames contain only 
small particles. That is, the large particles are most often found in concentrations 
so low that to capture more than one particle in a frame is rare. When there are 
more than five particles in an image frame, the local concentration is on the order 
of 100,000 L' 1 and higher. The small particles can be found at low concentration 
and at extremely high concentration. The conditional size distributions shown in 
Figure 19b could result from the aircraft flying first through a region of large 
particles in low concentration, followed by a region of small particles in high 
concentrations, similar to the large scale structure that caused the conditional 
spectra to separate in Figure 1 1 . However inspection of the CPI image frames 
themselves, shown in Figure 20, reveals that the regions of high concentrations 
of small particles are interspersed with regions of low concentrations of large 
particles, mostly bullet rosettes, on scales down to tens of meters and perhaps 
smaller. ' Baker et al. (2000) also observed regions with clumps of small particles 
separated by regions with single rosettes and aggregates of rosettes in a mid- 
latitude cirrus cloud. 

6. Summary 

This research investigated microphysical data collected by the NCAR C-130 
during the FI RE. ACE field experiment conducted over the Beaufort Sea in May 
and July 1998. The standard microphysical measurements in the NCAR C-130 
were supplemented, for the first time, with data collected by the cloud particle 
imager (CPI). CPI data were used to separate spherical images (i.e., water 
drops) from non-spherical images (ice particles) in mixed-phase clouds. The 
CPI images also provide a method for determining inhomogeneity in the cloud 
drop field. CPI data were combined with conventional PMS FSSP, 260X and 2D- 
C data to determine total particle concentration, and the phase discriminating 
capability of the CPI was used to determine the water and ice particle size 

distributions. 

A major focus of this investigation concentrated on the microphysical 
properties of Arctic boundary layer clouds, which were observed on 1 1 of the 
total of 16 missions flown by the C-130. Here, we define boundary layer clouds 
in the sense previously described by Curry et al. (1988), where the clouds are 
essentially low-lying stratus clouds that may or may not be thermodynamically 
connected with the surface. These boundary layer clouds were found to vary 
considerably, both from day to day and within the clouds themselves. The main 
microphysical features of the boundary layer clouds, based on data collected by 
the C-130, can summarized as follows: 

> From data collected during the 16 aircraft missions and the definition 
used here to define the presence of boundary layer cloud (i.e., thermally 
mixed to the surface and/or within 300 m of the surface), 11 of the 16 
days had boundary layer clouds. In May, six out of eight cases had 
boundary layer clouds and all six clouds were mixed from the surface to 




23 


cloud base. The depths of the mixing layers in May ranged from 1 50 to 
1200 m In July, five of the eight cases had boundary layer clouds, but 
none of these were mixed from the surface to cloud basa Thus, based 
on these (limited) data, low-lying clouds in the Arctic, which have been 
shown to strongly influence the melt of the Arctic Ocean ice pack (Curry 
et al. 1993), were prevalent and displayed significantly differe 
subcloud mixing characteristics in May and July. 

> Based on analysis of data from 21 vertical (slant) profiles flown i on four 

days in essentially all-water clouds, nine of the profiles m0 * J* 

homogeneous, adiabatic conditions with monomodal drop sze 
distributions and no drizzle drops. Data from the remaining 12 profiles 
indicate that these clouds were either too thin to have adequate LWC to 
“ne whether they were adiabatic, or that they were .actively 
mixinq from cloud top downward. In the latter cases, the LWC and 
temperature in the upper portion of cloud were generally non-adiaba ic 
and the clouds were inhomogeneous. The drop size distributions were 
often bimodal near cloud top and drizzle was sometimes observed. 

> The adiabatic regions of clouds were used to evaluate the " c ® 

of the LWC probes, although the adiabatic LWC never exceeded 0.4 g 
m~3 in these thin clouds, so the instruments were only evaluated within a 

, relatively low range of LWC. The two King 

exceed the adiabatic value and were within 75% of adiabatic l he 
Gerber PVM probe scattered around the adiabatic value and some im 
exceeded it by up to about 35%. The FSSP systematically exceeded 
the adiabatic LWC value by up to a factor of two. 

> A mixed-phase boundary layer cloud investigated on 4 May 1998 that 
was 360 m thick and ranged in temperature from -22 tc i 25 
displayed considerable variation in hydrometeor fields. Based on CPI 
measurements, regions about 10 Km across and separated by only 
100 m in the vertical contained either small cloud drops and few ce 
particles, drizzle, or graupel particles. Seeding from cirrus clouds aloft 
likely influenced the microphysical processes in this cloud. 

In addition to the study of boundary layer clouds, a deep stratus cl0 “ d 
cloud tese at 2 Km ( + 2 • C) and cloud top at 6 Km (-23 • C) was studied 
Examinattonof the microphysics in this deep stratus cloud revested extreme 
variability in cloud particles and two layers with exceptionally high (2,500 to 4,00 
L* 1 ) concentrations of small ice particles. 

Starting the description of the stratus cloud from the top down, a region of 

Sffl tmtJf-OT I", ice concentrations were observedf^coud 
top to the -13 0 C revel. In a thin layer around -12 ° C, very high (~ 2500 L ) ice 
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concentrations were observed, with about 97% of the ice identified as small (< ~ 
50 pm) spheroidal (but not perfectly round) particles, fragments and irregular 
shapes. Only 2% of the ice appeared to be vapor grown. Relatively low ice 
concentrations and small amounts of supercooled liquid water were observed 
between -1 1 0 C and -7 ° C. From -5.5 ° C to -3.3 ° C, another layer with very 
high (4000 L' 1 ) ice particle concentrations was observed. Although this layer was 
within the temperature regime where the Hallett-Mossop ice multiplication 
process is often observed, there were insufficient drops > 23 pm to meet the 
Hallett-Mossop criteria. Instead, the FSSP appeared to be responding to the 
small ice spheroids observed by the CPI, since the Rosemount icing detector did 
not measure any supercooled liquid water. The ice particles in this region were 
similar to those observed at -12 0 C, with the main difference being more 
columns and fewer crystals with side-plane growth. Small, spheroidal (but not 
perfectly round) ice particles, fragments and other irregular ice particles 
accounted for 85% of the total ice particle concentration. This is in general 
agreement with observations by Rangno and Hobbs (2000) in another FI RE. ACE 
stratus cloud, where they analyzed CPI images and found that 57% of the ice 
particles were “frozen drops" (presumably similar to our small spheroids), 
fragments and particles with irregular shapes. Rime breakup, fragmentation from 
crystal collision, and drop shattering are discussed as possible ice multiplication 
mechanisms, but there is no physical way to adequately verify any of these 
processes. Drizzle was observed precipitating through cloud base of this deep 
stratus cloud. 

An Arctic cirrus cloud was also investigated and shown to be highly 
inhomogeneous on scales down to tens of meters or less. Regions with high 
concentrations of small ice particles are interspersed with regions of low 
concentrations of large particles, mostly bullet rosettes. The CPI images 
showed that the regions with small ice particles were in very high (~1 00,000 L ) 
local concentrations, although the average concentration in the cirrus measured 
by the CPI and FSSP probes was on the order of several hundreds per liter. 

The variability in the microphysical properties of Arctic stratus (boundary 
layer deep stratus and cirrus) clouds, both within a cloud and from cloud to 
cloud, presents a challenge for Arctic column and process modelers alike. The 
variability of hydrometeor types and concentrations within Arctic stratus clouds 
also presents a challenge to investigators retrieving microphysical properties 
from remote measurements. While presenting these challenges, Arctic stratus 
clouds also offer the benefit of being relatively easy to study with large turboprop 
and jet aircraft in the summer months, in that they are persistent and cover 
extensive areas. 


Acknowledgments: We would like to thank Dr. Judith Curry for her insightful 
comments. We are also indebted to the NCAR Research Aviation Staff for their 
scientific and instrumentation support, and to the crew of the C-130 for their 




25 


skillful piloting and maintenance of the aircraft. This work has been funded under 
the NASA FIRE program, Contract No. NASI -9601 5. 


References 

Arnott P A. D. Mitchell, C. G. Schmitt, D. Kingsmill, and D. Ivanova Analysis of 
the FSSP performance for measurements of small crystal spectra in arms. 
Proceedings: 13 th International Conference on Clouds and Precipitation, 
Reno, NV., 14 - 18 August, 2000, 191 - 193, 2000. 


Baker B. A. R. P. Lawson, and C. G. Schmitt, Clumpy Cirrus. Proceedings. 
13 {h International Conference on Clouds and Precipitation, Reno, NV., 14 - 
August, 2000, 637 — 640, 2000. 

Baumgardner, D„ An analysis and comparison of five water droplet measuring 
instruments, J. Appl. Meteor., 22, B91-910, 1983. 

Baumoardner D W Strapp, and J. E. Dye, Evaluation of the forward scattering 
spedromeierVVobc: Part II. Corrections for coincidence and dead-t,me 
losses, J. Atmos. Oceanic Techno!., 2, 626-632, 1985. 

Baumgardner D„ Corrections for the response times of particle measuring 
probes, Sixth Symp. Meterol. Obs. and Instruments, New Orleans, Am. 
Meteorol. Soc., pp. 148-151, 1987. 

Baumoardner D and M. Spowart, Evaluation of the forward scattering 
spectrometer probe. Part Ilk Time r^ponse and laser "homogeneity 
limitations, J. Atmos. Oceanic Tech., 7, 666-672, 1990. 

Baumgardner, D„ W. A. Cooper, and J. E. Dye, Optol andekctronic jWtattons 
of the forward-scattering spectrometer probe Med ^ c 'g 90 
Measurements Techniques: 2nd Volume ASTM STP 1083, 115-127, 1990. 

Biter C J., J. E. Dye, D. Huffman and W.D. King, The drop-size response of me 
CSIRO liquid water content probe. J. Atmos. Oceanic Techno!., 4, 359-367, 

1987. 

Brenouier J-L Coincidence and dead-time corrections for particle counters. Part 
Twgh Concentration measurements with an FSSP, J. Atmos. Oceanic 
Techno!., 6, 585-598, 1989. 

Cerni T. A., Determination of the size and concentration of cloud drops with an 
FSSP, J ■ Clim. &Appl. Meteorol., 22, 1346-1355, 1983. 




26 


Cober, S.G., J.W. Strapp and G.A. Isaac, An example of supercooled drizzle 
drops formed through a collision-coalescence process. J. Appl. Meteor . , 35, 
2250-2260. 1996. 

Cober, S.G., G.A. Isaac, Korolev, A. V., and J. W. Strapp, Assessing the relative 
contributions of liquid and ice phases in winter clouds. Proceedings. 13 
International Conference on Clouds and Precipitation, Reno, NV„ 14 - 18 
August, 2000, 689 - 692, 2000. 

Cooper, W. A., Effects of coincidence on measurements with a forward scattering 
spectrometer probe, J. Oceanic Atmos. Tech., 5, 823-832, 1988. 

Curry, J.A., Interactions among turbulence, Radiation and microphysics in Arctic 
stratus clouds, J. Atmos. Sci ., 43, 90-106, 1986. 

Curry, J. A., E. E. Ebert and G. F. Herman, Mean and turbulence structure of the 
summertime Arctic cloudy boundary layer. Q. J. R. Meteorol. Soc. 114, 715- 
746, 1988. 

Curry, J.A., F.G. Meyer, L.F. Radke, C.A. Brock and E.E. Ebert, Occurrence and 
characteristics of lower tropospheric ice crystal in the Arctic, Int. J. Chmatol., 
10, 749-764, 1990. 

Curry, J.A., E.E. Ebert and J.L. Schramm, Impact of clouds on the surface 
radiation balance of the Arctic Ocean, Meteor. Atmos. Phys., 51, 197-217, 
1993. 

Curry, J.A., W.B. Rossow, D.A. Randall and J.L. Schramm, Overview of Arctic 
Cloud and Radiation Characteristics, J. of Climate, 9, 1731-1764, 1996. 

Curry, J.A., J O. Pinto, T. Benner and M. Tschudi, Evolution of the cloud 
boundary layer during the autumnal freezing of the Beaufort Sea, J. Geophys. 
Res., 102, 13851 - 13860, 1997. 

Curry, J.A., P. Hobbs, M. King, D. Randall, P. Minnis, et al., FIRE Arctic Clouds 
Experiment, Bull. Amer. Meteorol. Soc., 81, 5-29, 2000. 

Dye J.E. and D. Baumgardner, Evaluation of the forward scattering spectrometer 
probe. Part I. Electronic and optical studies, J. Atmos. Oceanic Technol. 1, 
329-344, 1984. 

Gardiner, B.A. and J. Hallett, Degradation of in-cloud forward scattering 
spectrometer probe measurements in the presence of ice crystals. J. Atmos. 
Ocean. Tech., 17T-180, 1985. 




27 


Garrett T. and P.V. Hobbs, Calibration of liquid water probes from the University 
of Washington’s CV-580 aircraft at the Canadian NRC wind tunne , Rep. 

Reseat Group, 20pp„ Dept, of Atmospheric Scences. 

Univ. of Washington, Seattle, WA, 1999. 

Gayet J G. Febvre and H. Larsen, The reliability of the PMS^SSP'nthe 
presence of small ice crystals. J. Atmos. Oceanic Techno!., 13, 1300-1310, 

1996. 

Gerber H B G Arends, and A.S. Ackerman, New microphysics sensor for 
airaaft use', Afmos. Res., 31, 235-252, 1994. 

Hallett, J., and S.C. Mossop, Production of secondary ice particles during the 
riming process. Nature, 249, 26-28. 1974. 

idpp climate* Chanae * The IPCC Scientific Assessment, edited by J.T. 
Houg“G ? jenkins and J.J. Ephramus, 365 pp„ Cambridge University 

Press, 1990. 

Jayaweera, K.O.L.F. and T. Ohtake, Concentration of ice crystals in Arctic stratus 
clouds, J. Rech. Atmos., 7, 199-207, 1973. 

Jensen J B P. H. Austin, M. B. Baker and A. M. Blyth, Turbulent mixing, 
spectral evolution and dynamics in a warm cumulus cloud. J. Atmos. Set., 42, 
173-192, 1985 

liana H G Feinqold W. R. Cotton and P. G. Duynkerke, Large-Eddy 
Qimi nations of Entrainment of Cloud Condensation Nuclei into the Arctic 
“aXenTMay^B FIRE7SHEBA Case Study, 2000 (This issue). 

King W D., D. A. Parkin and R J. Handsworth, A hot wire ■£« <*£ 
having fully calculable response charactenstics, J. Appl. Meteo ., , 

1813, 1978. 

Knollenberg, R.G., Techniques for 7’ Ho^fanTt 

Formation, Optical Properties, and Effects, edite y 
Deepak, Academic Press, New York, NY. 15-91, 1981. 

Korolev, A.V., J.W. Strapp and G.A. Isaac, Evaluation oMte Aaxiracyof PMS 
Optical Array Probes. J. Atmos. Oceanic Techno!., 15 , 708-720, 19 . 

Korolev AV G.A. Isaac and J. Hallett, Ice particle habits in Arctic clouds, 
Geophys.Res. Letters, 26, 9, 1299-1302, 1999. 

Laursen K., Project documentation summary, Surface Heat Budget of the Arctic 
Ocean, NCAR RAF Project No. 8-101, 1998 



28 


Lawson, R. P. and W. A. Cooper, Perf<^anre rf some airborne thermometers 
in clouds, J. Atmos. Oceanic Technol., 7, 480-494, 1990. 

Lawson R P and R H. Cormack, Theoretical design and preliminary tests of 

L two new particle spectrometers for cloud microphysics research. Atmos. Res., 
35, 315-348. 1995. 

Lawson R P Improved particle measurements in mixed phase and 

"canons on P climate modeling, Proceedings: WMO Wotohop on 
Measurement of Cloud Properties for Forecasts of Weather, Air Quality a 
Climate , 23-27 June, Mexico City, pp. 139-158, 199/. 

i axA/Qnn R P and A M Blyth A comparison of optical measurements of liquid 

La> water conten^and drop size distribution in adiabatic regions of Flonda cumuli, 
Atmos. Res., 47~48, 671-690, 1998. 

i r p and T L Jensen Improved microphysical observations in mixed 

phase clouds, Presented at 'the American Meteorological S°^ Con ^ ence 
on Cloud Physics, 17-22 August, -Everett, Washington, pp. 451^*54, 1998. 

I auucnn R P A V Korolev S G. Cober, T. Huang, J. W. Strapp and G.A. 

LaV Jsaac, improved Measurements of the Drop Size Distribution of a Freezing 
Drizzle Event, Atmos. Res., 47-48 , 181-191, 1998. 

Mo-rin i P A V Korolev A Heymsfield, G. A. Isaac, and S. G. Cober, 
Thermodynamics of icing cylinder for measurements^ liquid wa er content 
supercooled clouds. J. Atmos. Oceanic. Technol., 2000 (in press). 

Mossop, S. C„ Secondary ice particle production duringrime growth: The < effert 
of drop size distribution and rimer velocity, Q. J. R. Meteorol. S ., 

1124, 1985. 

MOSSOP, S. C„ and J. Hallett, J., Ice crystal coneentmtot ii n atmulus clouds: 
Influence of the drop spectrum, Science, 186 , 632-644, iy/«. 

.ar Jswjssaw ““ 

Year in ice gives climate insights, EOS, 80 (41), 484 , 

Pinto J. O., J. A. Curry, and J. M. Intrieri, Cloud-aerosol interactions dunng 
autumn over the Beaufort Sea, 2000 (This Issue). 



29 


Poellot M R W.P. Amott and J. Hallett, In Situ Observations of Contrail 
Microphysics and Implications for Their Radiative Impact. J. Geophys. Res., 
104 , (DIO), 12,077-12,084, 1999. 

Pruppacher, H. R. and J. D. Klett, Microphysics of Clouds and Precipitation, 
707pp., D. Reidel Publishing Company, 1978. 

Ranqno A L and P. V. Hobbs, Ice particles in stratiform clouds in the Arctic and 
possible mechanisms for the production of high ,ce concentrates, 2000 

(This issue). 

Rooers D. C. P. J. DeMott and S. M. Kreidenweis, Measurements of Ice 
Nucleating Aerosol Particles in the Arctic, 2000 (This issue). 

Vali, G., Ice multiplication by rime breakup. Communications i s Ma V^me 
conference intemationale sur la physique des nuages. Clermont-Ferrand, 
France, Volume I, 227-228, 1980. 

Wendisch, M., A, Keil, and A.V. Korolev, FSSP characterizati^ ^th 
monodisperse water droplets, J. Atmos. Oceanic Technol.. 13. 6, 1152 1 , 

1996. 

Wendisch, M. A, quantitative comparison of 9 ;ound^ FSSP and PVM 
measurements. J. Atmos. Oceanic Technol., 15, 887-900, 1998. 

Witte H J., Airborne observations of cloud particles and infrared flux density in 
the Arctic, MS Thesis, Dept, of Atmospheric Sciences, 102 pp., Umv. of 
Washington, Seattle, 1968. 

Yum Seong Soo, and James G. Hudson, Vertical distrtbutions of cloud 
condensation nuclei spectra over the spnngtime Arctic Ocean, 20 ( 

Issue). 

Zhang, CL, S. Krueger, Qiang Fu and Qing Wang, Large-Eddy Simulation of An 
Arctic Stable Cloudy Boundary Layer, 2000 (This issue). 



30 


Figure Captions 

FIGURE 1: Diagram showing the fundamental design of the CPI. See text for 
explanation of operation. 

altitude on the right. 

FIGURE 3: Histograms of the number -of Particles ;per fram 1 J^m'inFigure 2 

SESi*^ 

or more particles (black) and by using those particles in frames wi 
particles per frame (gray). 

FIGURE 4' A scatter plot of measured Liquid Water Contents (for four 
inSruments) taken during nine different slant profiles plotted against the 
calculated adiabatic values. 

FIGURE 5- The top panel shows the Liquid Water Content versus time and 

- " 

calculated adiabatic LWC profile assuming different cloud bases. 

FIGURE 6 The CPI 1WC and Concentration are shown along with the false 
LWC signal observed by the FSSP, PVM, and King probes. 

tin irf 7 Time series of the temperature and altitude recorded by the C-130 
^“sampled a 

adiabatic temperature are shown in P P series of LWC recorded by 

sl?e^ t0 the calcu,ated adiabatiC Va ' Ue in the 

bottom panel. 

Pin irf 8- Time series of measurements from a boundary layer cloud on 29 
ji?y U 1998 Conditional size spectra (as shown in figure 3) for each region a 

shown above. 

FIGURE 9- Time series of measurements and their standard deviations made by 
the C-1 30 as it passed through two boundary layer adiabatic clouds. 
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pin IRF 10- Time series of measurements and their standard deviations made 
byIhe^C-130 as it flew different level passes through boundary layer clouds on 
four different occasions. 

FIGURE 1 1 ' Histograms of the number of particles per frame for all the ■CPI 
imaged frames during the time period 22:09:00 - 22:09:40 shown in Figures 2 
and 9 9 along with the Poisson distribution with the same mean and (bottom) 
conditional 9 particle size distributions produced by using only those parties sin 
frames with P 5 or more particles (black) and by using those particles in fram 
with 3 or less particles per frame (gray). 

Fin irf 12 Time series of measurements and their standard deviations made 
by the C-130 Is it flew different level passes through a boundary layer cloud on 4 

May 1998. 

Pin irf IT ExamDle of a Boundary Layer Cloud with extremely variable 

hydrometeor fields that exist over relatively small (10 Km HE BA^ p 

, n flioht ^3^ j S for the NCAR C-130 as it descended over the SHEBA ship. 

Ibe ufs shl^ on 'the track correspond to the CPI derived 
distributions of ice and the FSSP size distributions of water ,n *^ e ™ x ® d _ ph 
clouds. Sample CPI images are shown for each partide size dis 

FIGURE 14- Concentrations measured by the FSSP 1 DC, and 2DC .are .plotted 
are shown to the right. 

FIGURE 15: Total particle size distributions measured by the FSSP O, CPI (soj'd 
> anri ?n C (M in the ‘central’ region with high ice concentrations near 12 C 
^rlm ediatlly above andbdow ^central r ? ,on. 

Water and ice particle size distributions are derived from CPI image data, 
examples of CPI images are also shown. 

FIGURE 16: As in Figure 15, except for the ‘central’ and ‘adjacent’ regions near - 
4.5 ° C. 

FIGURE 17: Particle size distributions measured by th ® FS ^Q^ ^rggp° lid 
line) and the 2D-C (A) for an arctic cirrus cloud observed on 29 J y 

FIGURE 18: Time series of particle concentrations measu^^ 

= SSSSW4KS5S5S SS 
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distribution with the same mean, and (bottom) conditional particle size 
distributions produced by using only those particles in frames wth 9 or more 
particles (black) and by using only particles that were imaged alone in a frame 

(gray). 

FIGURE 20: Images of arctic cirrus ice particles sampled on 29 July 1998 

durina FIRE ACE Each box with a large crystal and each group of small boxes 
represents one CPI image frame. The image frames reflect the relative positions 
oHhe particles along the flight track. A significant point here is that the spatial 
variation between low concentrations of large particles and high concentra ions 

of small particles 
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